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Field of the Invention 

[0003] The invention relates to the field of surface measurement and, more specifically, 

to the field of non-contact surface measurement. 

Background of the Invention 

[0004] Dimensional metrology, the measurement of the size and shape of objects, is very 

important in today's manufacturing environment in which machines perform much of the 
fabrication and assembly of complex objects composed of many subassemblies. The shape and 
size of each component in a complex assembly, such as an automobile, must be held to close 
tolerances to ensure that the components fit together properly. 

[0005] Ideally such measurements of shape and size are accomplished without physical 

contact in order to save time in making the measurement. Many non-contact measurement 



methods make use of available machine vision systems. The measurement of surface contour 
information is an especially difficult problem in machine vision systems since depth information 
is often lost or is difficult to interpret. To compensate for the loss of depth information and the 
difficulty in interpreting the information which is available, many machine vision systems utilize 
light to create moire pattems on the surface of the object in order to obtain contour information. 
One disadvantage of the moire technique is its inflexibility with respect to objects of different 
sizes. Objects of different sizes may require new corresponding physical setups. Due to this 
disadvantage, it is difficult to use a moire technique for large-scale objects. Another 
disadvantage of a moire technique is that the resolution obtained using the technique may not be 
high enough for many applications 

[0006] Literferometric methods have also been used when detailed measurements of the 

surface are needed. Although interferometric systems provide surface contour information, they 
are sensitive to vibrations in both the object being measured and the source of illumination being 
used. 

[0007] What is needed is a technique which avoids these problems by resolving fringe 

order ambiguities, mitigating degradation due to speckle effects and attaining high resolution. 

Summary of the Invention 

[0008] The invention is related to a method for determining, on an object having a 

surface, three-dimensional position information of a point on the surface of the object including 
the steps of providing two sources of radiation having a spatial distribution of spectral regions, 
illuminating the surface with radiation from each of the sources to produce a first fringe pattern 
at a first position on the surface, moving the first fringe pattern to a second position, generating a 
first wrapped cycle map in response to the first and second positions of the first fringe pattern, 
estimating fiinge nxunbers in the first fringe pattem, changing the first firinge pattem to generate a 
second fiinge pattem at a first position, moving the second fiinge pattem to a second position, 
generating a second wrapped cycle map in response to the first and second positions of the 
second fringe pattem, estimating fringe numbers in the second fringe pattem in response to the 
estimated fiinge numbers in the first fiinge pattem, and determining three-dimensional position 
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information of the point on the surface in response to the estimated fringe numbers in the second 
fringe pattern and the second wrapped cycle map. 

[0009] The invention is fiirther related to a method for synthesizing a wrapped cycle map 

corresponding to a fringe pattern of a predetermined spacing on a surface including the steps of 
providing two sources of radiation having a spatial distribution of spectral regions, illuminating 
the surface with radiation from each of the sources to produce a first fringe pattern at a first 
position on the surface, moving the first fringe pattem to a second position on the surface, 
generating a first wrapped cycle map in response to the first and second positions of the first 
fringe pattem, changing the first fiinge pattem to generate a second fringe pattem at a first 
position, moving the second fiinge pattem to a second position, generating a second wrapped 
cycle map in response to the first and second positions of the second fiinge pattem, subtracting 
the second wrapped cycle map from the first wrapped cycle map, and wrapping the difference 
between the second wrapped cycle map and the first wrapped cycle map to generate a wrapped 
cycle map corresponding to the fiinge pattem of the predetermined spacing. 

[0010] The invention fiirther relates to a method for determining, on an object having a 

surface, three-dimensional position information of a point on the surface of the object including 
the steps of generating a first and a second beam of radiation having a first frequency and being 
coherent with respect to one another, illuminating the surface with the second beam of radiation, 
generating a first interference pattem at a first position in response to the first beam of radiation 
and radiation from the second beam of radiation scattered by the surface, moving the first 
interference pattem to a second position, generating a first wrapped cycle map in response to the 
first and second positions of the first interference pattem, estimating intensity cycles in the first 
interference pattem, changing the first interference pattem to generate a second interference 
pattem at a first position, moving the second interference pattem to a second position, generating 
a second wrapped cycle map in response to the first and second positions of the second 
interference pattem, estimating intensity cycles in the second interference pattem in response to 
estimated intensity cycles in the first interference pattem, and calculating three-dimensional 
position information in response to the estimated intensity cycles in the second interference 
pattem and the second wrapped cycle map. 
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[0011] The invention further relates to a method for determining, on an object having a 

surface, three-dimensional position information of a point on the surface of the object including 
the steps of generating a first and a second beam of radiation being coherent with respect to one 
another, illuminating the surface with the second beam of radiation, generating a first 
interference pattem in response to the first beam of radiation and radiation from the second beam 
of radiation scattered by the surface, phase shifting one of the first and second beams of radiation 
with respect to the other of the first and second beams of radiation to generate a first phase 
shifted interference pattem, generating a first wrapped cycle map in response to the first 
interference pattem and the first phase shifted interference pattem, changing the first interference 
pattem to generate a second interference pattem, phase shifting one of the first and second beams 
of radiation with respect to the other of the first and second beams of radiation to generate a 
second phase shifted interference pattem, generating a second wrapped cycle map in response to 
the second interference pattem and the second phase shifted interference pattem, subtracting the 
second wrapped cycle map from the first wrapped cycle map, and wrapping the difference 
between the second wrapped cycle map and the first wrapped cycle map to generate a wrapped 
cycle map. 

[0012] The invention also relates to an apparatus for projecting fiinges onto a surface of 

an object including two sources of radiation having a spectral distribution, a collimator in optical 
communication with the two sources, the collimator generating two substantially coUimated 
beams of broadband radiation, a dififractive grating in optical communication with the collimator, 
and a lens in optical communication with the diffractive grating, the lens generating two images 
of radiation having a spatial distribution of spectral regions. 

[0013] The invention also relates to a method for projecting fringes onto a surface of an 

object including the steps of providing two sources of radiation separated by a distance and 
generating radiation having a spectral distribution, coUimating the radiation to generate two 
substantially coUimated beams of radiation, delineating spectral components of the coUimated 
beams of radiation, and generating two images of radiation having a spatial distribution of 
spectral components. 



-4- 



[0014] The invention further relates to a method for determining, on an object having a 

surface, three-dimensional position information of a point on the surface of the object including 
the steps of projecting a first cyclical pattern to produce a first projected pattem at a first position 
on the surface, moving the first projected pattem to a second position, generating a first wrapped 
cycle map in response to the first and second positions of the first projected pattem, estimating 
cycle numbers in the first projected pattem, changing the first projected pattem to generate a 
second projected pattem at a first position, moving the second projected pattem to a second 
position, generating a second wrapped cycle map in response to the first and second positions of 
the second projected pattem, estimating cycle numbers in the second projected pattem in 
response to the estimated cycle numbers in the first projected pattem, and determining the 
surface in response to the estimated cycle numbers in the second projected pattem and the second 
wrapped cycle map. 

[0015] The invention fiirther relates to a method for synthesizing a wrapped cycle map 

corresponding to a projected cyclical pattem of a predetermined spacing on a surface including 
the steps of projecting a first cyclical pattem to produce a first projected pattem at a first position 
on the surface, moving the first projected pattem to a second position on the surface, generating a 
first wrapped cycle map in response to the first and second positions of the first projected pattern, 
changing the first projected pattem to generate a second projected pattem at a first position, 
moving the second projected pattem to a second position, generating a second wrapped cycle 
map in response to the first and second positions of the second projected pattem, subtracting the 
second wrapped cycle map fi-om the first wrapped cycle map, and wrapping the difference 
between the second wrapped cycle map and the first wrapped cycle map to generate a wrapped 
cycle map corresponding to the projected cyclical pattem of the predetermined spacing. 

[0016] The invention fiirther relates to a method for mitigating the effects of speckle on a 

measurement of a point on a surface of an object including the steps of generating a coherent 
fiinge pattem, projecting the coherent fiinge pattem along an optical path onto the surface of the 
object such that the fiinge pattem substantially grazes the surface of the object, and detecting the 
fiinge pattem and the speckle in an image of the surface of the object, where a normal to the 
surface of the object is substantially orthogonal to the optical path. 
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[0017] The invention further relates to a method for mitigating the effects of speckle on a 

measurement of a point on a surface of an object including the steps of generating a coherent 
fringe pattern from two sources of radiation separated by a distance, projecting the coherent 
fringe pattern onto the surface of the object, detecting the fringe pattem and the speckle in an 
image of the surface of the object, translating the two sources such that the fringe pattem remains 
substantially stationary and the speckle changes, detecting the new fringe pattem and the changed 
speckle in an image of the surface of the object, and determining the fringe pattem substantially 
without speckle in response to the detected fringe patterns and the change in speckle. 

[0018] The invention also relates to a method for mitigating the effects of speckle on a 

measurement of a point on a surface of an object including the steps of generating a coherent 
fringe pattem from two sources of radiation separated by a distance, projecting the coherent 
fringe pattem onto the surface of the object, detecting the fringe pattem and the speckle in an 
image of the surface of the object, translating the object laterally with respect to the two sources 
of radiation along a path parallel to the equiphase plane of the fringe pattem such that the fringe 
pattem remains substantially stationary with respect to the surface of the object and the speckle 
changes, detecting the new fringe pattem and the changed speckle in an image of the surface of 
the object, and determining the fringe pattem substantially without speckle in response to the 
detected fringe patterns and the change in speckle. 

[0019] The invention also relates to a method for mitigating the effects of speckle on a 

measurement of a point on a surface of an object including the steps of projecting a coherent 
fringe pattem onto the surface of the object, providing a lens having a transmission fimction 
having a gradual transmission drop off at the lens edge, and detecting the fringe pattem in an 
image of the surface of the object, where the transmission frinction substantially reduces the 
speckle influence on the measurement. 

[0020] The invention also relates to a method for mitigating the effects of speckle on a 

measurement of a point on a surface of an object including the steps of generating a first fiinge 
pattem from two sources of radiation having a first frequency and separated by a distance, 
changing the first fiinge pattem to generate a second fiinge pattem by changing the first 
frequency to a second frequency, and changing the distance in response to the difference between 
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the first frequency and the second frequency, where the ratio of the distance to the difference 
between the first frequency and the second frequency is substantially constant. 



[0021] The invention further relates to a method for projecting fringes onto a surface of 

an object including the steps of providing two sources of radiation separated by a distance, 
having a spectral distribution, and being coherent with respect to the other of the sources, 
illuminating the point on the surface of the object with the radiation from each of the sources, 
moving one of the sources relative to the other of the sources, and detecting radiation scattered 
by the point on the surface of the object. In another embodiment of the invention, the step of 
providing two sources of radiation includes providing an initial beam of radiation having a 
spectral width, generating a first radiation beam at a first beam angle and a second radiation 
beam at a second beam angle from the initial beam of radiation, and imaging the first radiation 
beam and the second radiation beam to form the two sources of radiation. 

[0022] The invention also relates to an apparatus for projecting fiinges onto a surface of 

an object including two sources of radiation separated by a distance, each source having a 
spectral distribution, and being coherent with respect to the other of the sources, a control system 
moving each of the sources relative to the other of the sources, and a detector positioned to 
receive radiation scattered from the point on the surface of the object. In another embodiment of 
the invention, the two sources of radiation include, an initial source of a beam of radiation having 
a spectral width, a beam separator in optical communication with the initial source of a beam of 
radiation generating a first optical beam and a second optical beam, and an imaging system 
optically connected to the beam separator. The imaging system generates the two sources of 
radiation each respective source of radiation corresponding to the first optical beam and the 
second optical beam, respectively. 

[0023] The invention further relates to a method for determining, on an object having a 

surface, three-dimensional position information of a point on the surface of the object including 
the steps of providing two sources of radiation having a spectral distribution and being coherent 
with respect to the other of the two sources, providing a detector at the point on the surface, 
illuminating the point on the surface of the object with the radiation from each of the sources, 
moving each of the soiu-ces relative to each other, detecting the radiation at the point on the 
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surface of the object, and calculating position information in response to the movement of the 
sources and the radiation detected at the point on the surface of the object. 

[0024] The invention further relates to an apparatus for determining, on an object having 

a surface, three-dimensional position information of a point on the surface of the object including 
two sources of radiation having a spectral distribution and being coherent with respect to one 
another, a control system moving each of the sources relative to each other, a detector positioned 
at the point on the surface of the object to receive radiation illuminating the point on the surface 
of the object, and a processor receiving signals from the detector, the processor calculating 
position information of the point on the surface of the object in response to the movement of the 
sources and the radiation received at the point on the surface of the object. 

Brief Description of the Drawings 

[0025] This invention is pointed out with particularity in the appended claims. The above 

and further advantages of this invention may be better understood by referring to the following 
description taken in conjunction with the accompanying drawings, in which: 

[0026] Fig. 1 is a block diagram of an embodiment of the invention for making surface 

contour measurements; 

[0027] Fig. 2 is a block diagram of an embodiment of a system for producing the two 

sources of radiation shown in Fig. 1; 

[0028] Fig. 2a is a block diagram of another embodiment of a system for producing the 

two sources of radiation shown in Fig. 1; 

[0029] Fig. 2b is a block diagram of yet another embodiment of a system for producing 

the two sources of radiation shown in Fig. 1; 

[0030] Fig. 3 is a block diagram of an embodiment of apparatus for supporting the two 

sources of radiation of Fig. 1 at a fixed distance relative to one another; 

[0031] Fig. 4 is another embodiment of the imaging system in Fig. 1 . 
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[0032] Fig. 5 is a block diagram of an alternate embodiment of the invention for making 

surface contour measurements; 



[0033] Fig. 6 is a flow diagram of an embodiment of the steps utilized by the processor of 

Figs. 1 and 5 in making surface contour measurements; 

[0034] Fig. 6a is one embodiment of a portion of the flow diagram of Fig. 6; 

[0035] Fig. 6b is another embodiment of a portion of the flow diagram of Fig. 6; 

[0036] Fig. 6c is yet another embodiment of a portion of the flow diagram of Fig. 6. 

[0037] Fig. 7 is a block diagram of one embodiment of a detector and processor 

arrangement for use with the systems of Figs. 1 and 5; 

[0038] Fig. 7a is a block diagram of an alternate embodiment of a detector and processor 

arrangement including a multiprocessor for use with the systems of Figs. 1 and 5; 

[0039] Fig. 7b is a block diagram of another alternate embodiment of a detector and 

processor arrangement for use in the systems of Figs. 1 and 5; 

[0040] Fig. 8 is a block diagram of another embodiment of the invention for making 

surface contour measurements; 

[0041] Fig. 9 shows one embodiment of photodetector elements positioned on the surface 

of an object in accordance with the embodiment of Fig. 8; and 

[0042] Fig, 9a shows another embodiment of photodetector elements positioned on 

spring arms for use in the embodiment of Fig. 8. 

[0043] Fig. 10 is a block diagram of yet another embodiment of a system for producing 

the two sources of radiation shown in Fig. 1; 

[0044] Fig. 1 1 is a block diagram of yet another embodiment of a system for producing 

the two sources of radiation shown in Fig. 1; and 
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[0045] Figs. 12 and 12A are block diagrams of the embodiment of Fig. 1 . 

[0046] Fig. 13 illustrates an embodiment of a series of steps used to perform progressive 

fringe division. 

[0047] Fig. 14 is a flowchart of an embodiment of a series of steps used to perform 

progressive fringe division. 

[0048] Fig. 14a is a flowchart of an embodiment of steps used to perform progressive 

fringe division utilizing fringe synthesis. 

[0049] Fig. 1 5 is a flowchart of an embodiment of a method of progressive fringe 

division using fringe synthesis and frequency tuning. 

[0050] Fig. 16 illustrates a block diagram of an embodiment of the invention used in the 

method of Fig. 15. 

[0051] Fig. 17 illustrates an embodiment of a broadband interference-fringe projector of 

the present invention. 

[0052] Fig. 18 illustrates one embodiment of a liquid-crystal spatial light modulator for 

use with an embodiment of the invention. 

[0053] Fig. 19 illustrates one embodiment of two nested double slits for use with an 

embodiment of the invention. 

[0054] Fig. 20 illustrates an embodiment of the invention for generating accordion fringe 

motion. 

[0055] Fig. 21 illustrates an embodiment of a broadband interference-fringe projector of 

the present invention. 

[0056] Fig. 22 contains graphical representation of fringe number error AA^ as a function 

of ^/(2n) for an embodiment of the invention known as histogram-based derippling. 
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[0057] 



Fig. 22a is a graph of the slopes + 1 of the curves shown in Fig. 22. 
Description of the Preferred Embodiment 



[0058] While describing the embodiment of the invention, reference will be made to 

"sources" and "sources of radiation." These terms are meant to refer to any source of radiation, 
including highly localized sources of radiation. 

[0059] Referring to Fig. 1, and in brief overview, two sources of radiation Pi and P2 are 

separated by a fixed distance D and have spatial coordinates of (x\^\^\) and (^^21^2,^2), 
respectively. The radiation from each of the sources Pi and P2 is coherent with respect to the 
radiation from the other one of the sources. Each source, Pi and P2, directs its respective 
divergent beam of radiation 12 and 14 toward a point Po on the surface of an object 10. The 
distance from each respective source of radiation, Pi and P2, to the point on the surface Po is 
indicated by R\ and R2, respectively. i|; is the angle between the line extending from the origin to 
the point Po and the line extending between sources Pi and P25 6s is the angle between the z axis 
and the line extending between the sources Pi and P2, and a is the half angle subtended by the 
source points as viewed from Pq. Each beam 12, 14 is substantially polarized in the same 
direction as the other beam 14, 12 and may be independently scannable to simultaneously 
illuminate different regions on the object 10. Aheraatively, the entire object 10 may be 
illuminated simultaneously. 

[0060] Light scattered 20 by the point Po is detected by a photodetector 22. In one 

embodiment, the photodetector 22 comprises an array of photodetector elements providing a two 
dimensional image of the object 10 to be measured. In a further embodiment, the array of 
photodetector elements is a charge coupled device (CCD). The detector 22 provides an output 
signal 26 comprising one or more individual signals, each one associated with a corresponding 
one of the photodetector elements of the detector 22. 

[0061] In a preferred embodiment, a focusing element 24 is positioned between the point 

Po on the surface of the object 10, and the photodetector 22, so as to image the illuminated 
portion of the object including point Po onto the detector 22. Because of the roughness of the 
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surface of the object, and because the illuminating radiation is coherent, the focused image will 
be speckled. The output signal 26 from the photodetector 22 is the input signal to a processor 
unit 28. 

[0062] A polarizer 30, in one embodiment, is placed between the focusing element 24 

and the detector 22. Polarizer 30 is oriented in a direction to maximize its coincidence with the 
principal polarization component of the scattered light 20, so as to improve the speckle contrast 
or the contrast of the fringe pattem. With this arrangement, the signal-to-noise ratio associated 
with light scattered from the surface of the object 10 is maximized. 

[0063] In one embodiment, the processor 28 is a single processor which operates on 

detector output signals 26 associated with each of the photodetector elements of the detector 
array 22. In another embodiment, the processor 28 is a multiprocessor having a plurality of 
individual processors and each photodetector element provides an input signal to a respective one 
of the processors. In yet another embodiment, in which the detector 22 is a CCD array, a 
plurality of the CCD elements provide an input signal to a respective processor of a 
multiprocessor. With the multiprocessor arrangements, computations on signals from a plurality 
of individual photoelements occur substantially simultaneously, thereby enhancing the signal 
processing speed. 

[0064] A control unit 32 controls the operation of the sources of radiation. Pi and P2, so 

as to change the phase of the radiation of one of the sources relative to the phase of the radiation 
from the other source as measured at the point Pq on the surface of the object 10. The processor 
28 may be in communication with control unit 32 via signal line, or bus 34. For example, in 
certain applications it may be desirable for the processor 28 to process signals from the detector 
22 at specific times relative to the scanning of the sources Pi and P2 over the surface of the object 
10 or relative to the rate at which the frequency of the radiation from the sources is swept. Since 
such scanning and frequency sweeping operations are controlled by control unit 32, 
communication between the control unit 32 and the processor 28 is desirable in these 
circumstances. It will be appreciated that the control unit 32 and the processor 28 may be 
physically separate units or, altematively, may be implemented by a single processing system. 
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[0065] Referring now to Fig. 2, in one embodiment the sources of radiation Pi and P2 are 

formed from the radiation emitted from a tunable laser 40. The radiation beam 44 emitted by the 
tunable laser 40 is split by a beam splitter 48. The radiation beam 50 reflected by the beam 
splitter 48 is caused to diverge by a lens 52. The divergent beam is then reflected by a moveable 
aiming mirror 54. The radiation beam reflected by the aiming mirror 54 provides one of the 
sources of coherent radiation, Pi . Similarly, the radiation beam 46 passing through the beam 
splitter 48 is caused to diverge by a lens 58 which directs the divergent beam to a second 
moveable aiming mirror 60. The radiation beam reflected by mirror 60 provides the second 
source of radiation, P2. Aiming mirrors 54 and 62 may be pivotable to selectively illuminate the 
surface of object 10. They may also be moveable to vary the positions of sources Pi and P2. 

[0066] Referring to Fig. 2a, another embodiment of the sources of radiation Pi and P2 is 

shown to include a tunable laser source 40 providing a beam of radiation 44. The radiation beam 
44 passes through a lens 62 which causes the beam to diverge, providing divergent beam 64. 
Divergent beam 64 is then reflected by beam splitter 48 to provide a first beam 66. A second 
beam 68 passes through the beam splitter 48, as shown. Moveable aiming mirrors 54 and 60 
reflect beams 66 and 68 to provide sources Pi and P2, respectively. 

[0067] Referring to Fig. 2b, in another embodiment the sources of radiation Pi and P2 are 

split from the radiation emitted from a tunable laser 40 using a fiber optic splitter 56. Fibers may 
have beam-forming elements at their end to control or set the divergence angle of the two beams, 
and in one embodiment the beam-forming elements may be lenses. Sources Pi and P2 may 
alternatively be formed from a pair of tunable lasers which are frequency locked together. Other 
suitable embodiments of radiation sources include any sources which generate a wave having a 
controllable phase, such as microwaves and sonic waves. 

[0068] In one embodiment, the sources of radiation Pi and P2 are maintained at a fixed 

distance D from one another by attaching each source to one end of a bar comprised of a material 
having a small coefficient of expansion. In another embodiment, the sources of radiation Pi and 
P2 are not held at a fixed distance but instead the distance between them, D, is known to a high 
degree of accuracy. 
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[0069] One illustrative bar 70 for supporting radiation sources Pi and P2 at a fixed 

distance D relative to one another is shown in Fig. 3. A bar 70 is provided with sockets 74 at 
opposite ends thereof. A ball joint 76 is pivotally positioned within each of the sockets 74, as 
shown. Each of the ball joints 76 has an end of a fiber firom a fiber optic splitter 56 (shown in 
Fig. 2b) positioned therein and an aperture 80 through which divergent radiation passes. Fibers 
may have beam-forming elements at their end to control or set the divergence angle of the two 
beams and in one embodiment the beam forming elements are lenses. In operation, the ball 
joints 76 are pivotable as shown by arrows 78 within the respective socket 74 and may be under 
the control of control imit 32 (shown in Fig. 1). With this arrangement, the divergent beams 12 
and 14 provided by the sources Pi and P2 at the ends of the fibers can be directed as desired to 
illuminate all, or a portion, of the object 10 including the point Po to be processed, while 
maintaining a fixed separation distance D. 

[0070] Referring again to Fig. 1, the coordinates of point Po on the surface of object 10 

are (x,y^y Although the x and 3; coordinates of Po are generally directly determinable fi'om the 
geometry of the detector 22 and the object 10, taking into account any magnification by 
intervening focusing element 24, the depth coordinate z, where the z axis is defined as being 
parallel to the optical axis of the imaging system, is not directly obtainable. The depth 
coordinate, z, however can be measured by first considering the difference in path length 

5 = R2-Ri + So (1) 

fi-om the radiation sources Pi and P2 to the point Po on the surface of the object 10. The quantity 
So is included to account for any path length difference in the beams that may occur before they 
reach points Pi and P2. 

[0071] Ifs is non-zero, then changing the firequency of the radiation emitted fi:om sources 

Pi and P2 will result in the phase of the radiation firom one source, as measured at point Po, 
changing with respect to the other source. This phase change results in a modulation of intensity 
of the radiation at point Po. The change in firequency, Av, required to complete one cycle of a 
change in intensity is given by the expression: 
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Av = - (2) 



where c is the speed of light. Thus, by measuring the change in laser frequency, Av, needed to 
cause one oscillation of intensity, the path difference s may be determined. The measurement of 
z is then based on determining the value of s for each value of jc and;;, as discussed below. 

[0072] Improved accuracy in the determination of s is obtained by measuring Av over 

many oscillation cycles. In practice it is convenient to work in terms of the number of oscillation 
cycles N (not necessarily a whole number) induced by a total change in frequency B. 

[0073] N is given in terms of Av and B as : 

N = A (3) 
Av 

Elimination of Av from Eq. (3) using Eq. (2) yields the following expression for s in terms of N: 

s=-N (4) 
B 

[0074] It should be noted that the number of oscillation cycles N induced by a frequency 

scan of length B corresponds to the number of interference fringes that move across the 
measurement point Po during the frequency scan. The interference fringes are produced by the 
interference of radiation emitted from the sources Pi and P2. These oscillations occur whether or 
not the fringes are resolved by the imaging lens 24. 

[0075] An uncertainty AA^ in the measurement of N, corresponds to an uncertainty As in 3' 

of 

As^-AN^s— (5) 
B N 

[0076] Equation (5) indicates that if the uncertainty AJV to which a single oscillation cycle 

can be determined remains constant, the uncertainty As in 5 is reduced by a factor equal to the 
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number of cycles N that are measured. There are numerous methods for determining N to 
various levels of resolution AN that are known to those skilled in the art. Examples of methods 
yielding a resolution of roughly one oscillation-cycle count (AN=1) are to perform a fast Fourier 
transform (FFT) on the data sequence or to count zero crossings of the high-pass filtered signal. 
Improved resolution of a fi-action of an oscillation-cycle count (AN<1) can be achieved, for 
example, by finding the argument of the discrete Fourier transform (DFT) where the magnitude 
of the DFT is maximized or by inspecting the phase of the oscillation cycle at the ends of the 
fi-equency scan. One technique known to those skilled in the art for accurate inspection of the 
phase is to insert a phase modulator in one leg of the beam path, i.e., between the beam splitter or 
fiber-optic spUtter and one of the sources Pi or P2 in Figures 2, 2(a), and 2(b). 

[0077] If Ii, I2, and I3 are signal intensities corresponding to phase shifts induced by the 

phase modulator of -90°, 0*^, and 90°, respectively, then the phase <|) of the oscillation cycle is 
given by: 



An uncertainty A0 in phase is converted to an uncertainty AN in the cycle count by dividing by 

271. 



uncertainty of AN = 1 cycle, an uncertainty of As = 20 |am is provided. An uncertainty of AN = 
0.1 cycle would improve the xmcertainty in 5 to As = 2.0 jim, provided that the spread in s over 
the lateral resolution is smaller than this quantity. If the spread in s over the lateral resolution on 
the surface of the object is larger than As, then the improved resolution in the measurement of s 
may still result in an improved estimate of an average or representative value ofs over that lateral 
resolution. 

[0079] In terms of the coordinate system: 




(6) 



[0078] 



For a typical frequency scan of B=15THz for a tunable diode laser, and for an 



s = ^j(x-x2y +{y-y2Y +(2-22)' -^(x-^iY +{y-yiy +i^-hT (7) 
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To make the calculation simpler, assume that the two sources Pi and P2 are located 
symmetrically about the origin at (Xi^i^x) and i-x\,-yu-z\). Then Eq. (7) becomes, in terms of 




(8) 



Solving for z, Eq. (8) becomes: 



4(xc, +yy^)z^ ±^^l6ixx, +yy,y +(s' -4z^)is' -D^ -4x' -4y') 



(9) 



where D is the distance between the two sources Pi and P2. Thus z is determined to within an 

ambiguity due to the existence of the positive and negative roots of Eq. (9). One way to avoid 
this ambiguity is by illuminating the object 10 so that the ^=0 line (labeled 16 in Fig. 1 for the 
case 5o=0) does not bisect the region of the object to be imaged. One way of moving the s=0 line 
is to vary 5o in Eq. (1). 

[0080] The sensitivity of the system to changes in s is shown by the ratio of As/Az, where 

Az is the uncertainty in z introduced by an uncertainty As in the value of s. This ratio ranges 
between zero, for a system lacking any practical range sensitivity and two, for a theoretically 
maximal system. A value of two is impractical to achieve because the surface of the object 10 
would have to lie between the two point sources Pi and P2 and only one side of the surface could 
be illuminated from each beam. The ratio As/Az is calculated by taking the partial derivative ofs 
with respect to z, from which the following expression for the range resolution is obtained: 



In equation (10), Ro is the distance from the origin to Po and \|/ is the angle between the line 
extending from the origin to the point Po and the line extending from point Pi to point P2 as 



Az = As 




(10) 
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shown in Fig. 1 . A useful configuration that provides good range resolution is to set \|/=90° for 
which the expression for Az simpHfies to 



Az= ^ (11) 

2 sin a cos 

where 9s and a are as shown in Fig. 1. In terms of Ro and D, tana=D/(2Ro). Equation (1 1) 

shows that the range resolution improves as the angle a increases and the angle 0s decreases. For 
values of As=5|am, a=10°, and 9s=45°, the range resolution is Az=20|im. 

[0081] Uncertainties (Ax, Ay) in the lateral position (x, y) of point Pq also affect the range 

resolution Az. If the two source points lie in the x-z plane, then the measurement of z is 
insensitive to uncertainties Ay, For v|/=90°, uncertainties Ax in jc cause an uncertainty 

Az = Axtan^, (12) 

in the measurement of z. Therefore, angles near 0s=O*^ offer the best immunity to uncertainty in 
the lateral position of point Pq. 

[0082] Because the depth of focus decreases as the lateral resolution of the optical system 

improves, there is a tradeoff between lateral resolution and maximum object depth. One method 
for reducing this limitation in object depth is to sequentially focus on different range planes and 
use only those pixels that are within the depth of focus. For example, a 100)Lim lateral resolution 
would limit the depth of field to the order of 1 cm, and an object with a 10 cm range could be 
imaged at full resolution by focusing sequentially at ten different ranges. To minimize the effects 
of depth of field, the z axis can be defined in a direction that minimizes the range extent, i.e., 
normal to the average plane of the surface of the object. To increase the lateral area that can be 
imaged without losing lateral resolution, multiple cameras (i.e., detector arrays 22) can be used to 
cover the entire area of interest of the object 10 or individual cameras can be used for inspecting 
regions of interest. Alternatively, the focal plane of single lenses can be populated with a 
plurality of detector arrays. These arrays can be translated independently to inspect various 
regions of the object at high resolution. Translation of individual detector arrays along the z axis 
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or tilting of the detector arrays can achieve simultaneous focusing for regions of the object at 
different depths to increase the allowable object depth. 



[0083] A potential difficulty with the optical imaging system in Fig. 1 is that the bistatic 

angle between the sources and the detector may introduce shadowing effects. These effects can 

be reduced by placing the lens closer to the sources as in Fig. 4 and using the lens in an off-axis 
configuration where the detector is offset laterally in the image plane. If the lens is designed for 
this purpose or has a sufficiently large field of view, then aberrations resulting firom off-axis 
imaging can be minimized. 

[0084] Referring to Fig. 5, an alternate embodiment of the present invention includes a 

moveable radiation source Pi and a stationary radiation source P2, each providing a divergent 
beam 150 and 154 and having a path length labeled R\ and R2 between such radiation source and 
a point Po on the surface of an object 10, respectively. The sources Pi and P2 may be generated 
by any suitable source of coherent radiation, such as a monochromatic laser, which is split to 
provide the two point sources Pi and P2. Moreover, various techniques are suitable for spUtting 
the radiation firom the coherent radiation soxu-ce, such as the beam splitter embodiments of 
Figs. 2 and 2a and the fiber-optic splitter embodiment of Fig. 2b. 

[0085] The divergent beams 150 and 154 are directed toward a surface of an object 10 on 

which a point Po is located having position information which is to be measured. Illxmiination 
scattered by the surface of the object 10 is focused by a focusing element, or lens 158 to impinge 
on a detector array 22. The lens can be used in an off-axis configuration as illustrated in Fig. 4 to 
reduce shadowing effects due to the bistatic angle. An optional polarizer (not shown) of the type 
described above in conjunction with Fig. 1 may be positioned between the focusing element 158 
and the detector array 22 in order to improve the contrast of the speckle image or the fringe 
pattern incident on the detector array 22. 

[0086] The detector array 22 is in communication with a processor xmit 28 for processing 

the image incident on the detector, as will be described. A control unit 32 is in communication 
with at least the moveable source Pi for moving the source Pi along an axis 160. As noted 
above, the control unit 32 and the processor unit 28 may be implemented by separate devices or 
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alternatively, may be part of a single system. Additionally, the control unit 32 and the processor 
unit 28 may communicate with each other, as may be desirable in certain applications. 



[0087] As described above in conjunction with Fig. 1, the depth coordinate z associated 

with a point Po on the surface of the object 10 can be determined as a function of the difference, 
R2-R1, between the path lengths Ri and R2 of beams 150 and 154, from sources Pi and P2 
respectively, to point Pq. hi the embodiment of Fig. 5, the phase of the radiation from moveable 
source Pi is changed by moving the source Pi along the axis 160 under the control of the control 
unit 32. With this arrangement, oscillations in the intensity at point Po are produced. 



[0088] The instantaneous coordinates of moveable point source Pi are 



x\-al^^ y\- ants, and z\ = aris 



(13) 



where a represents the magnitude of translation of point source Pi. and ^s? and ris are 
direction cosines representing the direction of translation with respect to the x, y and z axes, 
respectively. The phase difference of the radiation from the sources Pi and P2 as measured after 
propagation to point Po is given by: 



(14) 



where (|)o represents a constant phase offset that may exist between the two coherent sources Pi 
and P2. As Pi translates along axis 160, the value of Ri changes, causing if to vary as a ftmction 
of a. 



[0089] The number of intensity oscillations that occur at point Po (or interference fringes 

that cross point Po) as source Pi moves away from the origin is given by: 



Ik X X 



yjx' +/ +z' -4{x-at,y+iy-am^y +(z-a«J^ 



(15) 



where Ro is the distance between point Po and the origin of the coordinate system, (|)(a) is the 
optical phase difference in Eq. (14) for a source separation of a, and (|)(0) is the optical phase 
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difference in Eq. (14) for a equal to 0. Consideration of Eq. (15) reveals that the number of 
intensity oscillations, N, resulting from movement from source Pi is independent of the location 
of the stationary source P2. This independence permits the sources Pi and P2 to be positioned in 
close proximity to one another. With this arrangement, the divergent beams 150 and 154 from 
respective sources Pi and P2 experience common disturbances, such as air turbulence and 
vibrations. In this way, the effects of such disturbances are minimized. Additionally, beams 150 
and 154 reach the surface of the object 10 with substantially identical polarization. 

[0090] Note also that a small separation between the sources Pi and P2 produces a fringe 

pattern on the surface of the object 10 with a large fringe spacing that can be easily resolved by 
the imaging lens 24. 

[0091] Since the magnitude of translation a of point source Pi is relatively small as 

compared to the value of Ro, Eq. (15) can be approximated to second order in a/Ro as follows: 



a 

N=- 



1.2^1 . 2 

cos^ — sm^ cos^^sm y^— 

2 K ^ ^0 J 



(16) 



where \\f is the angle between the line extending from the origin to the point Pq and the line 
defined by the direction of translation of Pi. 

[0092] Eq. (16) indicates that to lowest order in u/Rq, knowledge of N allows the angle \\f 

to be determined. Given knowledge of v|/ from three or more locations, the (x^^) coordinates of 
Po could be determined through triangulation. We now describe an embodiment similar to the 
one corresponding to Fig. 1, where the x mdy coordinates are determined from location of the 
image point in the detector array. 

[0093] The measurement of z for a given (x,y) location can be made either by determining 

the number of intensity oscillation cycles N corresponding to a movement of Pi by a distance or 
by measuring the rate at which such intensity oscillations occur. Consider first a measurement of 
2 based on determining the number of cycles N. With N known, all of the variables in Eq. (15) 
are known except for z. Solving Eq. (15) for z yields the following expression: 
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_^_ An^.pp'-(r-n',)(x'^/) ^^^^ 



where 



A = xi,+ym,+^(p'-l) (18) 



and 



-P-^ (.9) 
a 

Equation (19) defines a dimensionless parameter having a magnitude ranging between zero and 
unity that represents the average modulation rate of the speckle intensity in terms of oscillation 
cycles N per wavelength unit traveled by Pi. For values of a approaching zero, Eq. (17) can be 
approximated as: 



z = 



—2 2 

P -n. 



[0094] The expressions for z in Eqs. 17 and 20 can be simplified by setting = 0, so that 

the translation of source Pi is confined to the x-y plane. This arrangement represents a good 
practical choice for translation of source Pi, as described below. The resulting expression for z 
can be written as follows: 

z= V^o (21) 

where the distance Ro from the scattering point Po to the origin of the x, y coordinate system is 
given by the exact expression: 

Ro = zr^ (22) 
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When a is small, Ro can be approximated as: 



Ro=^^^4^ (23) 
P 

[0095] Consider now the measurement of z based on knowledge of the instantaneous rate 

at which the intensity oscillations occur. The instantaneous oscillation rate p can be expressed in 
a manner similar to the average oscillation rate in Eq. (19), as follows: 

p = X— (24) 
aa 

Substituting the expression for the number of intensity oscillations, N, from Eq. (15) into 
Eq. (24) yields: 

P= I 2 2 ^^^^ 

^j(x-a£^y '^(y-'am^Y +(z-a«,)^ 
where the relation: 

ls^ + ^ws^ = «s^ = l (26) 
has been used to simpUfy the numerator. For small values of a, p can be approximated as: 



p«cos^-sm I// cos^sm^ — ^ (27) 

Rq 2 Rq 



Solving, Eq. (25) for z yields: 

^^—2 T (28) 

P -< 

where 

0 = [x£^^ym^+a(p'-l)]n^ (29) 
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and 



When Ws = 0, Eq. (29) can be written in the form of Eq. (21), with: 



(30) 



_ ^a[p^ -l)[2(^^, ^ymya\^(xl^ ^ym^f 

Ro n v31) 

\P\ 

For small values of a, Eqs. (28) and (31) can be approximated by Eqs. (20) and (23), 
respectively, with p replaced by p. 

[0096] In order to estimate range resolution, consider the imcertainty Az in the 

measurement of z that would be introduced by an uncertainty (AN or Ap) in the quantity being 
measured. For simplicity, this calculation is based on the approximate expression for N given by 
Eq. (16). To find Az, we take the partial derivative of N (or p) with respect to z and equate this 
derivative to the ratio AN/Az (or Ap/Az), to yield: 

Az = G^JLAN = GRoAp (32) 
a 

where 

G= 7 i ^ = . ^ . ry (33) 

""(^^5 +nn^)n sin^lcos^^sin^-cos^sin^^cos(^-^^)l 

is a geometrical factor that accounts for the direction of translation and the direction to the 
scattering point, hi the first form for G, 

^=x/Ro, m=y/RQy and«=z/Ro (34) 

are direction cosines for the Point Pq. hi the second form for G, 0 and ^ are the polar and 
azimuthal angles, respectively, representing the direction firom the origin to Pq in a spherical- 
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polar coordinate system. Likewise, the direction of translation of the source point is given by 0s 
and (|)s. 

[0097] Consideration of Eq. (32) reveals that range resolution degrades with increasing 

object distance Ro and improves with increasing magnitude of translation a of source Pi. 
Consideration of Eq. (33) reveals that the geometrical factor G ranges between unity and infinity, 
where unity corresponds to the best range resolution achievable. 

[0098] The optimal direction of translation of source Pi for a given scattering-point 

direction is obtained firom Eq. (33) by choosing £ s, ^s, and «s such that G is minimized for the 
given values ofi^m and «. Application of this constraint yields: 

(^,,m,,/zj= =-_=,Vl-n' (35) 

^ Vl-«^ Vl-«^ J 

which implies that the optimal translation direction is orthogonal to the line extending from the 
origin to the scattering point Po (v|/=90°) and lies in the plane of incidence formed by said line and 
the z axis ((l)s=<|)). Substitution of the values in Eq. (35) into Eq. (33) results in: 

G = -^=L== ^ =^ (36) 
Vl-«' V^'+m' sm^ 

From Eq, (36), it is observed that the best achievable G value of unity occurs when n = 0 {Q = 
90°), which implies that the scattering point lies in the x-y plane. It is also observed that the 
resolution degrades such that G approaches infinity for scattering points lying on the z axis. For 
example, G = 2 for 9 = 30° and G = 5.76 for 9 = 10°. Although it is not possible to satisfy 
Eq. (35) for every point in the image without changing the translation direction for each point, 
the condition for optimal resolution can be approximated by satisfying Eq. (35) for a 
representative image point. 

[0099] By Eqs. (25) and (27), the instantaneous modulation rate p depends on the offset 

magnitude a of the translating point. For techniques based on measuring p, it is desirable for p to 
vary as little as possible during the scan so that there is nearly a one-to-one correspondence 
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between values of /? and z. Then standard spectral-analysis techniques can be applied to estimate 
the value of p and determine 2. To quantify the degree of nonuniformity in p that occurs during a 
scan, we define: 

^ ^ p{ci) - piO) ^ pja) - cos y/ ^^^^ 
p(0) cosy/ 

Substitution of the approximate form for p from Eq. (27) into Eq. (37) and keeping only the 
lowest order term containing a, yields: 

(38) 

cost^ 

Equation (38) states that the modulation nonuniformity increases linearly in the ratio a/Ro of scan 
length to object distance. Furthermore, the nonuniformity vanishes when v|/=0° and increases 
without bound when \|/=90®. We observe, however, that there is no range resolution when \|/=0° 
because all points on the \|/=0® line have the same modulation rate, regardless of range, i.e., G=oo 
in Eq. (33). Therefore, there is a tradeoff between minimizing the nonuniformity and obtaining 
optimal range resolution. 

[0100] A reasonable measurement configuration that simultaneously provides good range 

resolution and reduced modulation nonuniformity is to set ns=0 and to use an off-axis optical 
system with the offset in the direction, i.e., <|)=<j>s. Then Eq. (33) for G reduces to: 

G= ~^ (39) 
sin(2^) 

[0101] As an illustrative example of the measurement technique, suppose it is desired to 

image an object that is 200 mm by 200 mm in the x-y plane fi'om a distance of Ro=lm using a 
laser of wavelength A.=0.7|am. If ns=0 and the center of the object is located at 0=30° and ^=^s, 
then, by Eq. (39), the geometric factor G will vary between 2.1 and 2.6 over the field of view. By 
Eq. (32), a translation of a=5nim will produce a range uncertainty of Az=16)im (in the middle of 
the image) for an uncertainty in the number of oscillations of one-twentieth of a count, i.e.. 
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AN=0.05. The total number of oscillation counts for the entire scan is N=3600 by Eq. (16). To 
estimate the modulation nonuniformity at the center of the image, we set \|/=60*^ in Eq. (38) and 
obtain x=0.0075 so that there is less than a 1% nonuniformity over the scan. This nonuniformity 
could be reduced further by introducing slight variations in the scan rate during the scan to 
compensate for any change in frequency during the measurement. 

[0102] Fig. 6 depicts an illustrative series of steps to be executed by the processor 28 of 

Figs. 1 and 5 to determine the depth coordinate z at each point (jc, y) on the object. The processor 
begins (step 100) by measuring a parameter of the intensity of the radiation scattered by a 
plurality of illuminated points on the object surface (step 108). From this information, the z 
coordinate for each measured point is calculated (step 112). 

[0103] An optional filtering process may be performed in step 116. Suitable filters 

known to those skilled in the art include, but are not limited to, smoothing filters, median filters 
and curve-fitting filters. Thereafter, the mapped points can be displayed or output in anyway 
known to one of ordinary skill in the art, following which the process is terminated in step 124, 
as shown. In one embodiment, the mapped points are plotted as a function of the computed z 
information on a mesh plot in step 120. 

[0104] Referring also to Fig. 6a, one embodiment of steps 108 and 112, suitable for use 

with the embodiment of Fig. 1 is shown. In step 108', the intensity of the scattered illumination 
is measured as a fiinction of laser fi-equency offset and N is measured using one of the methods 
known to those skilled in the art. Thereafter, s is calculated for each location (x,y) in step 110' 
using Eq. (4) and z is calculated for each location (x,y) in step 112* using Eq. (9). 

[0105] An altemate embodiment of process steps 108 and 1 12 for use in conjunction with 

the embodiment of Fig. 5 is shown in Fig. 6b. In this case, the parameter of the intensity 
measured in step 108" is the number of times, N (not necessarily a whole number), that the 
intensity cycles as the moveable source Pi (Fig. 5) translates. Once N has been determined in 
step 108" through one of the methods known to be skilled in the art it is converted to p by 
Eq. (19) in step 110'*. z is then calculated in step 1 12" with the use Eqs. (17) and (18). Another 
embodiment of process steps 108 and 1 12 for use in conjxmction with the embodiment of Fig. 5 
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is shown in Fig. 6C. Here, the parameter of the intensity measured in step 108" is the 
instantaneous oscillation rate p at which oscillations occur as the source point Pi translates, p is 
converted to z in step 1 12'" through Eqs. (28)-(30). 

[0106] Various arrangements of detector 22 and processor 28 are possible. In one 

embodiment, shown in Fig. 7, the photodetector elements 22i,i to 22n,m of the detector array 22 
are read out serially. The serial output 36 of detector array 22 provides an input to a processor 
28. Processor 28 may comprise a single processor or alternatively, may be a multiprocessor 
comprising a plurality of processors. 

[0107] Referring also to Fig. 7a, an alternate detector and processor arrangement is 

shown. In this embodiment, the processor 28 is a multiprocessor comprising a plurality of 
processors 28i,i to 28n,nT Each of the photodetector elements 22ij to 22n,m of the detector array 
22 provides a respective output signal 38 to a corresponding one of the processors 28i,i to 28n,m. 
With this arrangement, each of the processors 28i,i to 28n,m is able to operate substantially 
simultaneously, in order to provide substantial performance advantages. More particularly, each 
processor 28i,i to 28n,m in the multiprocessor unit 28 is responsible for making the z coordinate 
calculation based upon the data received from the corresponding element 22i,i to 22n,m of the 
photodetector array 22. Thus, the z coordinate for each location of the surface of the object 10 
may be determined rapidly. 

[0108] Fig. 7b shows a further altemate embodiment of the detector and processor 

components for use in the systems of Figs. 1 and 5 in the form of a unitary detector and processor 
array 25. The array 25 is fabricated on, and supported by, a common substrate or is fabricated as 
a Multi-Chip Module (MCM) or with Surface Mount Technology (SMT). The detector portion 
of the array 25 includes photodetector elements 22i,i to 22n,m and the multiprocessor portion of 
the array includes processors 28i,i to 28n,m- More particularly, each of the detectors 22i,i to 22n,ni 
is associated with, and positioned adjacent to, a respective one of the processors 28i,i to 28n,m 
and provides an input signal to the respective processor, as shown. The processors 28i,i to 28n,m 
process the information from the respective detectors 22i,i to 22n,m substantially simultaneously 
to provide the determined depth coordinates. 
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[0109] Referring to Fig. 8, another embodiment of the invention includes an array of 

detectors 22' placed against the surface of the object 10 whose surface contour is to be measured. 
With this arrangement, rather than observing the light scattered from the point Pq on the surface 
of the object 10 to determine z, the measurement of the phase shifting of the light is performed 
directly at the surface of the object. Although not shown, the system of Fig. 8 includes a control 
unit 28 for controlling sources Pi and P2 and a processor 28 for processing radiation incident on 
detector 22' as shown and described above in conjunction with Figs. 1, 5 and 6. 

[0110] The arrangement and mechanism for locating the photodetector elements 23 on 

the surface of the object 10 may vary. In one embodiment shown in Fig. 9, a plurality of 
individual photodetector elements 23 of array 22' are positioned on the surface of the object 10 in 
the area of interest. 

[0111] In another embodiment, shown in Fig. 9a, individual photodetector elements 23 of 

the array 22' are mounted on spring arms 84 cantilevered from a support and control unit 88. The 
spring arms 84 are moved over the surface of the object 10 by the control unit 88 in order to 
contact specific points, or regions of interest. The cantilevered support of spring arms 84 causes 
each individual detector 23 to remain in contact with a location on the surface of the object 10 as 
the arms 84 are moved thereover. That is, as the contour of the object surface varies, the spring 
arms 84 move up and down accordingly. 

[0112] It will be appreciated that two or more additional radiation sources may be used in 

apparatus and methods of the present invention. For example, an additional source or sources 
can be used to determine x,y coordinate information regarding the object or a portion thereof 
Additionally, extra radiation sources may be used to reduce any processing inaccuracies or 
ambiguities attributable to shadowing of a region of interest. 

[0113] It will also be appreciated that other variations to the embodiment involving 

moving source points may be used. For example, the two points may both move with opposing 
motion, they may both move in the same direction with constant separation, they may rotate 
about a common center point, or motion may be simulated by using an array of source points that 
can be switched on and off by the control system. 
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[01 14] Movement of both sources can provide advantages over a single moveable source 

such as reduction of the modulation nonuniformity in Eq. (38). A strategy for reducing this 
nonuniformity is to eliminate the linear term in Eq. (27) for r by moving the points in such a way 
that N is an odd function of a. Splitting the translation equally between the two sources so that 
their center remains fixed (referred to as opposing motion) results in the following expression: 



a 



\2 



X — /, 



\2 



(40) 



The expansion of N corresponding to Bq.(16) that contains terms up to fourth order in a is: 



X 



1 . 2 a 

cosT — cosTsin ^— r 
8 R' 



(41) 



Thus, the second-order term in a and all higher-order even terms in a are eliminated and the 

magnitude of the third-order term is decreased by a factor of four. 



[0115] 



An exact expression for z based on Eq. (17) and Eq. (22) is thus given by: 



K +yf"sK ±pJk +yf"sf -[p^ -«.^) /4)|p^ -ij 

Z = : 



-2 2 
P -"s 



(42) 



When Us = 0, Eq. (42) reduces to Eq. (21) with 



(43) 



For small values of a, Eqs. (42) and (43) reduce to Eqs. (20) and (23), respectively. Elimination 
of the first-order dependence on a by moving both sources, however, makes this approximation 
more accurate than for a single moving source. 



-30- 



[0116] For opposing motion of two sources, the expression for the modulation rate given 

in Eq. (25) becomes: 




2 a y f a y f a 
x + — L\ + v + — /wJ + z +— 
K 2 V V 2 V V 2 ' 



Using a small a approximation that is valid up to and including third order of a, Eq. (44) 
becomes: 



3 

p « cos4^ — cos^sin' 4^— (45) 

8 Rn 



In comparison with Eq. (27), the absence of the linear term in Eq. (45) results in an average 
oscillation rate p that remains more constant during a scan. The corresponding modulation 
nonimiformity % previously defined in Eq. (37) can now be expressed as 

jc« — sin^^— (46) 

Eq. (46) yields a much smaller value for x than Eq. (38) because of the second-order dependence 
on the ratio a/Ro and because of the absence of cos \|/ in the denominator. One of the advantages 
of moving both source points is that the nonuniformity remains small for the optimal range- 
resolution situation where \|/=90°. 

[01 1 7] The results for range uncertainty Az previously given for motion of a single source 

are based on the first term in the small-a approximations of N and p. Because the first terms in 



-31- 



these approximations are identical for motion of two sources, the previously described results are 
also valid for the motion of two sources. 



[01 18] Referring to Fig. 10, another embodiment of the present invention includes two 

oppositely moveable sources Pi and P2. A beam 44 from a laser source 40 is split into a 
transmitted beam 46 and a reflected beam 50 by beam splitter 48. The transmitted beam 46 is 
reflected from a mirror 94 to a lens 58. Similarly, the reflected beam 50 is reflected from a 
mirror 96 to a lens 52. A right angle prism 92 with mirrored surfaces 91 and 93 is used to 
redirect the beams 46 and 50 from the lenses 58 and 52, respectively. The beams 46 and 50 
converge at points Pi and P2, respectively, and continue as divergent beams 82 and 86, 
respectively. Phase shifting can be accomplished by changing the optical path length difference 
between the beam paths followed by the separated beams 46 and 50. This can be achieved, for 
example, by inducing out-of-plane translations of at least one of the elements 48, 94, or 96 with a 
piezoelectric transducer (not shown). Altemately, a phase-shifting device (not shown) may be 
inserted in one of the beam paths 46 or 50. An example of such a device is a rotating glass disk 
that has been coated such that its optical thickness varies in discrete steps with rotation angle. 
Prism 92 is translatable along an axis 99 to affect a symmetric motion of the sources Pi and P2. 
Alternatively, the prism 92 can remain fixed while the laser 40, beamsplitter 48, mirrors 94 and 
96 and lenses 52 and 58 are translated as one group along the axis 99 in order to cause a 
symmetric translation of the sources Pi and P2. Changing the separation of the sources Pi and P2 
can also be accomplished using an angular modulation device (not shown) positioned in the 
beam 44 between the laser 40 and beam splitter 48 which causes a change in the beam angle. 
Examples of angular modulation devices are galvanometer mirrors and acousto-optic modulators. 

[0119] hi another embodiment (not shown), the sources Pi and P2 are virtual sources. 

That is, the sources Pi and P2 are not small sources of light or regions where the beams 46 and 50 
pass through focus. Instead, diverging beams 82 and 86 are generated using optical components 
which cause the beams 46 and 50 to diverge (i.e. negative lenses). The resuhing illumination at 
the object can be indistinguishable from other embodiments using real sources Pi and P2. 

[0120] Fig. 1 1 shows yet another configuration for generating oppositely moving sources 

Pi and P2. The beam 44 from a laser source 40 is split into a transmitted beam 46 and a reflected 
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beam 50 by beam splitter cube 48. The transmitted beam 46 is reflected from a curved mirror 95 
and from a mirrored surface 91 of a right angle prism 92. Similarly, the reflected beam 50 is 
reflected from a curved mirror 97 and from a mirrored surface 93 of the right angle prism 92. 
The beam splitter cube 48 can be modified (e.g., removal of glass outside the beam path) to 
permit close mounting of optical components. The curved mirrors 95,97 perform the functions 
of the flat mirrors 94,96 and lenses 58,52 in Fig. 10. The curved mirrors 95,97 can be off-axis 
parabolic mirror segments which yield small focal spots at Pi and P2. Alternatively, the curved 
mirrors 95,97 can be spherical elements. Aberrations caused by spherical elements will not 
affect the fringe contrast, but only the fiinge position. These aberrations can be incorporated into 
the analysis. Phase shifting can be accomplished by translating beam splitter 48 or by placing a 
phase-shifting device (not shown) in one of the beam paths 46 or 50 as discussed previously. 

[0121] Referring again for a moment to Fig. 1, the derivation in the prior description is 

based on a telecentric imaging configuration. The lateral position of the image point Pi in the 
focal plane is related to the lateral position of the object point Pq by the magnification of the lens 
24 and is independent of the separation between the lens 24 and the object 10. Defocusing blurs 
the image Pj but does not affect its apparent lateral position at the detector 22. Consequently, the 
lateral position of an object point Pq can be estimated by the position of its image point Pi, 
regardless of the distance. Unfortunately, telecentric imagers 24 may be impractical for large 
objects because the lens 24 must be at least as large in diameter as the largest dimension of the 
measured object 10. Thus most lenses are not telecentric. 

[0122] Referring now to Fig. 12, to avoid this problem, a more general approach for 

determining the (x^y^) coordinate of the point Po is to assume that the position of the image point 
Pi at the detector 22 is a measure of the angle ax between the object point Pq and the optical axis 
98 of the focusing element 24. The coordinate can be derived from the fiinge number N and the 
location of the image point at the detector 22. The following equations are exact (do not assume 
small a). The lateral coordinates are given by: 

X = jc^ + (2 - z,)tanaj^ (47) 

and 
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y = yi+{z-z,)tanaj 



(48) 



where and are the components of angle a as measured in the x-z and y-z planes, 
respectively. 



[0123] 



The z coordinate is given by: 



p'b-b,c,±a 



z = z«. + - 




(49) 



where 



A = 1 + tan^ ttjj + tan^ Uy 



(50) 
(51) 



(52) 



Bs=ljtana^ \anay 



(53) 
(54) 
(55) 



and 



3'p=r/-:>'m-^/tana. 



(56) 



The coordinate (Xm, ym, Zm) represents the mid-position Pm between points Pi and P2 and allows 
for a generalized solution with (xm, ym, Zm) located at other than the origin of the coordinate 
system. For increased versatility, Pm and other parameters may vary with time. 
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[0124] The generalized determination of object coordinates described above can be 

implemented in software on a computer which receives data from a processor unit in 
communication with a plurality of detectors 22 (e.g., a CCD array). Alternatively, the algorithm 
can be implemented directly in the processor unit 28. 

[0125] The use of Eqs. (47)-(49) to determine the (x, y, z) coordinates of a point Po 

requires knowledge of the angles ax and ay in addition to the normalized fringe number ~p as 
defined by Eq. (19). In the example above, these angles are determined by the location of the 
image point P,- in the focal plane 22. In order to determine ax and ay accurately, it is necessary to 
correct for any image distortion that may occur in the lens 24. Distortion can be corrected, for 
example, by creating a distortion function that maps angles ax and ay onto image points Pi. The 
distortion map may vary for different lenses and may even vary slightly for the same lens 24 in 
different imaging configurations. 

[0126] In another embodiment of the present invention, the inconvenience of 

characterizing lens distortion is avoided by using additional source heads to determine the 
remaining two direction components necessary for triangulation. For example, as depicted in 
Fig. 12a, an additional source head is placed at the position Pi and the lens 24 which was 
previously at Pi (see Fig. 12) is placed at any convenient position for viewing the region of 
interest on object 10. If the sources Pix and P2x in the second source head are oriented parallel to 
the X axis, then a measurement of the fiinge number Nx at point Po for this source head is 
equivalent to measuring the angle v|/x between the lines defined by points Pix and Pi, and points Pi 
and Po (see Eq. (41)). The angle ax required for Eqs. (47)-(49) is then determined from \}/x by the 
relationship ax = 7r/2 - \|/x. Likewise, (not shown in Fig. 12a) the angle ay can be determined by 
using a source head at position Pi with sources Piy and P2y oriented parallel to the y axis. The 
fringe number Ny for this orientation of a source head provides a measurement of the angle ay 
through ay = 7c/2 - \|/y. The source heads providing the angles ax and ay can, for example, be the 
same source head rotated 90° about the z axis or different source heads made to appear at the 
same location through use of a beam splitter. 
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[0127] When triangulating using multiple source heads, it is not necessary to know the 

location of lens 24 and detector 22 (or the location components of detectors 23 placed on the 
object 10). The purpose of the imaging system 22 and 24 (or detectors 23) is to observe the 
illumination intensity at Pq. Each pixel in the detector array 22 represents an observation point 
Pq. Consequently, an (x, y, z) coordinate may be determined for each pixel. Various approaches 
may be used to separate the measurements of N corresponding to different source heads. For 
example, the measurements may be time multiplexed by illuminating sequentially with different 
source heads or wavelength multiplexed by simultaneously illuminating with multiple source 
heads operating at different wavelengths. Measurements can be further distinguished through the 
use of polarization diversity. 

[0128] To further elucidate the process of triangulation using multiple source heads, 

reference is made again to Eq. (41). For typical apphcations, a is small compared with Ro and 
Eq. (41) can be simplified to. 



7\^ = -cos^ = — COST (57) 
A c 

[0129] The dependence on source spacing a and wavelength X, or frequency v, is 

eliminated from Eq. (57) by writing N in terms of the normalized fringe number p using 
Eq. (19): 

p = cos4^ (58) 

Thus, the normalized fringe number corresponding to a particular source head is a direct measure 
of a component \\f of direction from the mid-point Pm of the source head to the observed point Pq. 
The three-dimensional coordinates of Pq can be determined through triangulation given 
measurements from three or more source heads 

[0130] Suppose that source heads are located at the three positions Pmi, Pm2, and Pm3, and 

the orientations of the source pairs corresponding to each of these locations are given by the 
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direction cosines (Iji, nisi, n si), (ls2, ins2, Hs:), and (ls3, niss, nss). In terms of Cartesian 
coordinates and direction cosines, Eq. (58) takes the form 



Pi = 



A^-^'n)' ^[y-y^n)' -^[^-^.S 



(59) 



where the subscript / distinguishes between source heads. Determination of Po coordinates (x, y, 
z) is achieved by measuring the normaUzed fringe numbers p i, P2 and, for each source head 
and solving the set of three simultaneous equations represented by Eq. (59). The formalism of 
Eq. (59) is convenient because it handles generalized source-head placement and orientation in a 
compact form. 

[0131] For those cases where the object point Po corresponding to a particular pixel 22 is 

not illuminated by three or more source heads, the missing direction information for that pixel 22 
can be filled in by interpolating between the directions calculated for surroxmding pixels 22. 
Furthermore, direction information obtained firom multiple source heads can be combined with 
direction information obtained from pixel location to improve measurement quality. 

[0132] The fact that multiple source heads makes it possible to determine Po without 

knowing receiver parameters (such as position, orientation, magnification, and distortion) makes 
it a simple matter to combine information from multiple imagers 22 and 24 (only one shown for 
clarity). This fact follows because the point clouds produced for each camera are already in the 
same coordinate system and do not need to be registered by going through a calibration 
procedure. Multiple imagers 22 and 24 can be positioned to view specific regions of interest at 
any desired view angle and lateral resolution. Thus, one imager 22 and 24 can give an overview 
of the entire object 10, and other imagers 22 and 24 can inspect regions of interest at high 



Po on the surface of the object 10 is determined by counting intensity-oscillations caused by 
motion of one or both of the sources Pi and P2. In its most rudimentary form, cycle counting 



resolution. 



[0133] 



In one of the above major classes of embodiments, the fiinge number N at a point 
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requires sampling the intensity at least twice per cycle in order to keep track of the count. Also, 
to minimize cost, it may be desirable to count cycles using standard off-the-shelf components 
such as a charge-coupled device (CCD) connected to a frame grabber. Because the number of 
counts increases the farther a point Po in the image frame is from the N = 0 position, it may be 
necessary to acquire and process hundreds of image frames using this approach. 

[0134] The need for acquiring many image frames to determine N can be overcome by 

using a more sophisticated counting technique. This technique, denoted as progressive fringe 
division, takes advantage of the periodic and predictable nature of the intensity modulation. 
Thus, N can be determined by sampling the intensity modulation at a few discrete values of the 
source spacing a, 

[0135] As depicted in Fig. 13, progressive fringe division consists of varying the fringe 

size from an initial coarse spacing (Col. 1) to a final fine spacing (Col. 3) in discrete steps. At 
each step, a wrapped phase map, or cycle map, (Row 4) is generated from phase-shifted fringe 
patterns (Rows 1-3). Each successive wrapped phase map is unwrapped based on the previous 
unwrapped phase map in a bootstrapping operation. The ratio of fringe spacings between steps is 
typically of the order of 10 or more, allowing a very dense fringe pattem to be unwrapped using a 
small number of data frames. Calculating the position coordinates from a dense fringe pattem 
improves precision. One of the advantages of progressive fiinge division is that the algorithm 
correctly assigns the fringe order number to each fringe, regardless of the complexity of the 
surface of object 10. For example, object 10 can contain surface discontinuities, large cliffs, and 
holes, or consist of a number of disjoint objects. Because progressive fiinge division does not 
rely on surface continuity, it can be used to locate the position of a number of unconnected 
discrete points that are distributed in space. 

[0136] A ftirther explanation of progressive fiinge division and its embodiment is aided 

by considering three facts. (For simplicity, opposing motion of the sources Pi and P2 is assumed 
in this discussion, although the following results can be formulated for other situations, including 
motion of only one source Pi and stationary frequency-tuned sources.) First, the number of 
intensity oscillations N that occur at a point Po as the source spacing increases from an initial 
value of zero to the final value a is equal to the fiinge number N at Po for the final value a. For 
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example, if 10 were a simple object such as a plane, one could count N fringes on the surface of 
the object between the N = 0 position (occurring at Ri = R2) and the point Pq. The value of N 
includes partial cycles or fringes and thus is not necessarily a whole number. 

[0137] This first fact can be further explained by noting that N can also be interpreted as 

the number of fringes that move past the point Po to be collected inside the region between Po and 
N = 0 as the source spacing a increases from zero to its final value. In other words, when a-0 
there are no fringes within this region, and the number of fringes contained inside this region 
increases by one every time a fringe moves across Pq. This correspondence between fringe 
number N and oscillation cycles is key to the robust three-dimensional imaging technique and 
makes it possible to determine N at an arbitrary point Po without making assumptions about the 
continuity of the surface surrounding the point Pq. Consequently, discontinuities in the surface 
such as holes and cliffs (i.e., edges) are handled without ambiguity. 

[0138] A second fact is that N is directly proportional to the separation a according to 

Eq. 57. Consequently, it is not necessary to complete an entire scan with a running from zero to 
its final value to determine N(a). The value of N for a spacing a2 can be estimated from its value 
at a spacing a\ using the relationship 

Nestimate(«2) = (^2/^1) N(fli) (60) 

Likewise, 

Nestimate(^JE2--fll) = N(a2) -N(ai) (61) 

which states that knowledge of the value of N at two source spacings a\ and a2 can be used to 
estimate its value at the smaller difference spacing a2-a\. 

[0139] A third fact is that a phase-shifting measurement can be made at any intermediate 

value of the source spacing a to provide a high-resolution map of the partial-cycle or * Vrapped" 
component of the fringe number N at each measurement point Pq within the field of view. In 
doing so, the wrapped component is calculated from the phase by 
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Nwrapped = ^/(2Tt) 



(62) 



and the phase is determined, for example, by Eq. (6). Thus, if (j) ranges between 

±71, Nwrapped ranges between ±0.5. Although a phase-shifting measurement can determine Nwrapped 

to high precision, it leaves the integer- valued offset for N undetermined. 

[0140] Because finer fiinges represent a smaller variation in p per fringe, they can 

provide higher precision in the determination of p (assuming successful unwrapping). For 
example, finer fiinges reduce the effect that phase-step errors have on the measurement of the 
surface profile. Finer fringes, however, also make it more difficult to perform the unwrapping 
because of the increased fiinge density. The advantage of progressive fi-inge division is that it 
allows very fine fiinges to be unwrapped correctly to achieve the highest possible resolution. 

[0141] To implement progressive fiinge division, wrapped cycle maps (as depicted in 

Row 4 of Fig. 13) are produced for the imaged surface at each of a number of discrete source 
separations a. It is often sufficient to produce these cycle maps for three or fewer fiinge sizes: 
coarse, medium, and fine. Successive unwrapped cycle maps are used in a bootstrapping 
operation to unwrap the next cycle map without introducing errors in the integer offset. The 
unwrapped cycle map corresponding to the finest fiinge pattem provides the highest-resolution 
measurement. 

[0142] Fig. 14 depicts an illustrative series of steps used to execute progressive fiinge 

division. The procedure begins (step 150) by setting the source spacing at an initial value a \ 
(step 154). This value is normally the smallest in a series of source spacings and is chosen so as 
to produce a coarse fiinge pattem on the object 10. A wrapped cycle map is then produced (step 
158) by performing a phase-stepping measurement, as illustrated by Eq. (6). Step 158 is 
completed by converting phase to cycles by Eq. (62). 

[0143] Independently of step 158, (step 162) an initial unwrapped estimate of N for the 

coarse fiinge pattem corresponding to ^ i is obtained. Step 162 may be performed 
simultaneously with step 158, before step 158, or after step 158. However, both step 158 and 
162 must be performed in order to continue to the next step. There are many approaches for 
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producing this initial estimate. For example, a\ can be made to be so small that the entire object 
falls within one fringe (e.g., -0.5 < N < 0.5) this making unwrapping unnecessary. Alternately, a 
can be large enough to produce a small number of fringes across the object 10 so that a simple 
model (such as a linear ramp) for the variation of N over the field of view will suffice as an 
estimate for N as long as this estimate is accurate to within one-half of a cycle at all measurement 
points to be unwrapped. Another example of a method for obtaining the initial estimate is to 
assume a generic shape or an expected shape and to calculate N based on Eq. (40). Other 
examples are to use the result of a previous measurement on the same or a similar object 10, or to 
use cycle counting to partition the field of view into cycles. The number of frames necessary for 
cycle coxmting is small here because of the small number of fiinges corresponding to ai. 

[0144] The purpose of the next step (step 166) is to unwrap the wrapped cycle map 

obtained from the measurement in step 158. If the unwrapped estimate from step 162 is accurate 
to within one-half of a cycle, unwrapping is achieved by use of the expression 

N(a) = Nwrapped(a) + R0Und[Nestimate(a) - Nwrapped(a)] (63) 

[0145] Note that the second of the two terms on the right-hand side of Eq. (63) is simply 

the integer offset required for unwrapping. Subtracting the wrapped measiurement from the 
unwrapped estimate in this term produces a stair-step-like ftmction whose value is approximately 
equal to the integer offset. Rounding to the nearest integer makes the expression exact, as long 
as the error falls within the one-half-cycle limits. Otherwise Eq. (63) introduces an integer offset 
error in N that may vary across the field of view. 

[0146] The result of step 166 is an improved estimate of N for the smallest source 

spacing a\. In order to fiirther improve the quality of the measurement, a new source spacing a2 
is set in step 170 that produces a finer fiinge pattem on the object 10. A new phase-stepping 
measurement is then completed in step 174 to yield a new wrapped cycle map corresponding to 
the finer fiinge pattem. The unwrapping of this measurement requires a new unwrapped estimate 
corresponding to a2. This estimate is obtained in step 178 by using Eq. (60) to scale the previous 
unwrapped estimate (obtained for the source spacing a\ from step 166) into the appropriate 
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range. Once scaling has taken place, the new wrapped cycle map from step 174 is unwrapped in 
step 182 by substituting it and the scaled estimate from step 178 into Eq. (63). 

[0147] The resuh of the first pass through steps 170 through 182 is a refined estimate of 

N that corresponds to the larger source spacing ai. Steps 170 through 182 can now be repeated 
with successively finer fiinges (larger source spacings) until optimal p resolution is achieved. 

(Each time through the loop the source-spacing subscript is incremented by one.) Normally, 
optimal resolution is achieved in one to three iterations. Typically, the ratio between the new and 
old source spacings is of the order of ten or more. The maximum value this ratio can take is 
determined in practice by requiring that the error between the estimated and true value of N 
remains under the one-half cycle limit. 

The loop represented by steps 170 through 182 in Fig. 14 is exited once the desired number of 
iterations has been completed. In step 186, N is converted to ~p through Eq. (19). If the 
objective of the measurement is to determine fiiU three-dimensional coordinates, a triangulation 
takes place using, for example, Eqs. (47)-(49) or Eq. (59), to obtain coordinates for each 
measurement point (step 190). Step 194 represents the end of progressive fiinge division. It 
should be appreciated that Fig. 14 is only an illustration of the procedure to be followed for 
progressive fiinge division, and that the component steps can be executed in a number of 
different orders. For example, it may be desirable to collect all data before beginning the 
calculations. 

[0148] The ability to produce fringes of widely varying size for use with progressive 

fiinge division requires the source spacing to vary by the same ratio as the fiinge spacing. A 
technique is now described for synthesizing a fiinge pattem of a desired spacing from two other 
measurements. One advantage of this approach is that the two measurements used to synthesize 
the fiinge pattem may be more readily obtainable. 

[0149] Fringe synthesis is based on the following modification of Eq. (61): 

Nwrapped(^2-^l) = [Nwrapped(«2) -Nwrapped(fll)] wrapped (64) 
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This result follows from wrapping both sides of Eq. (61) and observing that wrapping the 
difference N(a2) -N(ai) on the right-hand side of Eq. (61) is equivalent to wrapping the 
difference Nwrapped(«2) -Nwrapped(«i) of individually wrapped components. In other words the 
difference between the wrapped quantities within the square brackets in Eq. (64) ranges between 
±1 and this difference must be wrapped to He within the correct range of ±0.5. 

[0150] Equation (64) provides a means for synthesizing phase-stepping measurements of 

wrapped N at a source spacing equal to the difference az-ai from individual phase-stepping 
measurements obtained at source spacings a\ and 02- This equation has a remarkable property for 
small values of a2-a\. Although the individual measurements from which the synthesized 
measurement is composed may be highly irregular and contain numerous discontinuities, the left- 
hand side of Eq. (64) can be smooth over the field of view, having only the number of wrapping 
discontinuities that correspond to the difference a2-ci\. Thus fringe synthesis is an effective 
technique for producing an initial phase map for beginning progressive fringe division. Fringe 
synthesis can also be used, to synthesize any of the remaining unwrapped phase maps in the 
sequence to be used for progressive fringe division. If p is to be calculated from a synthesized 
cycle map, Eq. (19) is modified to read 

P = —~-^ (65) 
Aa 

where Aa represents the spacing difference between the two measurements. 

[0151] Fig. 14a illustrates one way that fringe synthesis can be incorporated into the 

progressive-fringe-division flow chart in Fig. 14. Step 154' is a modification of step 154 to 
include the selection of a base separation ao in addition to the separation ai. Step 158'. is a 
modification of step 158 to measure wrapped cycle maps at both separation values ao and ai. 
Steps 160' (inserted after Step 158) and 176' (inserted after step 174) are additional steps that 
make use of Eq. (64) to synthesize the wrapped cycle map corresponding to the differences an- 
ao. Step 186' is a modification of step 186 that allows p to be calculated from a synthesized 
cycle map using Eq. (65). 
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[0152] One of the advantages of fringe synthesis is that it may reduce the demands on a 

system for producing motion between sources Pi and P2. For example, a small incremental 
change from the widest source separation can produce fringes corresponding to a very small 
source separation. This approach not only reduces the travel requirement, but also overcomes 
some potential limitations of the soiu-ce heads described in Figs. 10 and 11. For example, with 
these source heads it is difficult to produce extremely small source separations, and fiinges may 
become unstable due to vibrations when attempting to do so. Another advantage of fringe 
synthesis is that it may improve measurement repeatability and accuracy due to the fact that it is 
normally easier to repeat a constant difference between two source separations than it is to 
maintain an absolute source-spacing calibration. 

[0153] Yet another advantage of fringe synthesis is that it minimizes errors of motion due 

to source-head misalignment. These motion errors may occur, for example, if the focal spots Pi 
and P2 in Figs. 10 and 1 1 do not overlap for the nominal a = 0 position. There are two types of 
misalignment to be considered. First, suppose that one of the focal spots Pi is displaced out of 
the plane of the diagram with respect to the other P2 by an amoimt d. Even though these two 
focal spots Pi and P2 each follow a straight-line trajectory and these trajectories can be parallel, 
the two spots miss by the distance d at their closest approach (i.e., the nominal a = 0 position). 
As the spots Pi and P2 pass at closest approach, the orientation of the fringe pattern on the object 
10 flips by 90° with respect to the nominal, or desired, orientation. Furthermore, the fringe 
spacing is never larger than that corresponding to a source separation d. (Alignment of the 
source head is facilitated by adjusting the displacement d until the fiinges stay properly oriented 
near a =0.) 

[0154] In the second type of misalignment, one of the focal spots Pi is displaced along 

the direction of motion 99 with respect to the other focal spot P2. This type of misalignment is 
caused by a sideways displacement of prism 92 (i.e., the displacement is perpendicular to axis 99 
and is in the plane of Figs. 10 and 1 1). At closest approach, one source Pi Hes directly behind the 
other source P2, producing a circularly concentric fiinge pattern on the object. As a increases, the 
fringes exhibit less curvature and approach their normal appearance. (Minimization of this 
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source of error during source-head alignment is facilitated by adjusting the prism 92 sideways for 
small a to minimize fringe curvature.) 



[0155] The difiFiculty that motion errors create for progressive fringe division is due to the 

resulting phase errors which are magnified when Eq. (60) is used to produce a scaled estimate of 
N. Motion errors limit the magnitude of the a ratio that can be used without introducing half- 
cycle errors in the estimate of N. Fringe synthesis overcomes this limitation because motion- 
induced phase errors are common to both measurements and cancel in the subtraction. Thus, 
even though the measurements may produce fringes that are distorted from the desired shape, the 
synthesized fringes are normal. 

[0156] Cancellation of phase errors between the two beams with fringe synthesis has the 

additional advantage of automatically compensating for a global phase-offset error between the 
two beams. These phase-offset errors are introduced, for example, by a path-length difference So 
between the two separated beams that is not an integer multiple of the wavelength X. Phase offset 
errors cause the N = 0 position to occur at a fringe phase of other than 0° (i.e., at other than an 
intensity maximum). These phase errors automatically cancel when performing fringe synthesis, 
making it uimecessary to calibrate the phase of the fringe pattem. Altematively, phase 
calibration can be achieved by determining the precise N = 0 position by fringe synthesis and 
adjusting the phase-shifting element for constructive interference at that point. 

[0157] The cancellation of phase errors for synthesized fringes is important for the 

additional reason that it tends to cancel out wavefront errors in the illumination pattems falling 
on object 10 from the individual sources Pi and P2. Wavefront errors arise from the use of 
imperfect or extended sources and produce perturbations or ripples in the fringe pattem. These 
perturbations translate into small height errors in the measurement. Cancellation of wavefront 
errors occurs in fringe synthesis because these errors are largely independent of the distance 
between sources. Therefore, the error is common to the measurements taken at each source 
spacing and cancels due to the subtraction in Eq. (64). Fringe synthesis is an effective means for 
reducing errors due to illumination defects and can serve to decrease system cost by reducing 
requirements on the quality of the source-head optics. 
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[0158] Fringe synthesis allows progressive fringe division to be carried out with a 

reduced range of motion of the source points thus allowing a simpler source head design. Recall 
that in general, a source head provides two types of fringe modulation — a variation in fringe size 
and sideways motion of the fringe pattern arising from phase shifting. For a pure phase shift, the 
phase is changed without altering the position of the source. Another way of changing phase, 
however, is to move one source Pi towards or away from the object relative to the other source 
P2. Although doing so changes the location of the source Pi, this change is very small (less than 
a wavelength), and the resulting change in orientation of the source pair Pi and P2 is small if the 
sources are sufficiently far apart. Because the source separation can remain large for fiinge 
synthesis, both types of modulation can be achieved through small motions of one or both source 
points Pi and P2. This reduced motion requirement allows simple source heads to be designed. 
For example, a source head might consist of the fiber-optic-splitter arrangement depicted in 
Fig. 2b. Phase shifting and variation in fiinge size are accomplished by moving the end of one of 
the fibers in orthogonal directions. 

[0159] Fringe synthesis has an important additional advantage of enabling progressive 

fiinge division to be used for measurements based on laser- frequency tuning. Without fiinge 
synthesis, application of progressive fiinge division to the measurements would require an 
extremely high fractional tuning bandwidth in order to achieve the large variations in fiinge size 
that are necessary for progressive fiinge division. Fringe synthesis makes it possible to produce 
fiinge-spacing variations that Ue far outside the conventional tuning range. The timing approach 
can provide measurements with very high repeatability and precision because of the fact that 
source motion is eliminated and frequency can be measured and controlled to high precision. 

[0160] Application of progressive fiinge division and fiinge synthesis with frequency- 

tuning is explained with respect to Eq. (57). N has the same fiinctional dependence on laser 
frequency v as it does on source spacing a, (It is assumed that the path-length difference 
represented by So in Eq. (1) is zero, otherwise this offset must be accoimted for in the equations.) 
Therefore, Eqs. (60), (63), and (64), which form the basis for progressive fiinge division and 
fiinge synthesis, also apply to frequency tuning if N is written as a fijnction of frequency v. The 
expressions that govem frequency tuning are then 
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Nestimate(V2) = (Va/Vi) N(Vi) 



(66) 



Nwrapped(V2-Vi) = [Nwrapped(V2) -Nwrapped(Vl)]wrapped (67) 

and 

N(V) = Nwrapped(v) + Round[Nestimate(v) - Nwrapped(v)] (68) 

Equations (19) and (65) also take the modified forms 



-p-^^ (69) 

a V 



and 



^ = £i^ (70) 

a Av 

respectively, where Av represents the difference frequency if N has been synthesized from two 
measurements. 

[0161] Fig. 15 illustrates the steps to be followed for execution of progressive fringe 

division using fringe synthesis and frequency tuning. (Note that Fig. 15 also appUes to 
progressive fringe division using fringe synthesis and source motion if v is replaced by a,) In this 
illustration, all data is collected at the beginning of the process. The procedure begins (step 200) 
by selecting the laser frequencies vo, vi, . . . Vn (step 204) to be used in the measurement of 
wrapped cycle maps Nwrapped(vo), Nwrapped(vi), . . . Nwrapped(Vn) (step 208). Equation (67) is then 
used to synthesize wrapped cycle maps for the difference frequencies vi-vo, V2-vo, .... Vn-vo 
(step 212). An estimate of the unwrapped cycle map is then obtained (step 216) for the first 
difference Vi-vo using one of many possible approaches (e.g., step 162 of Fig. 14). A preferred 
approach is to begin with a difference Vi-vo that is small enough so that N(vi-vo) lies within the 
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±0.5-cycle limit. In step 220, Eq. (68) is used to calculate the unwrapped cycle map for the next 
larger difference Vi-vo in the series. This unwrapped cycle map is scaled using Eq. (66) in step 
224 to be fed back into step 220 to unwrap the next synthesized cycle map in the series. The 
loop is exited after the cycle map corresponding to the final frequency difference Vn-vo has been 
unwrapped. In step 228, Eq. (70) is used to convert the final unwrapped cycle map N to fiinge 
nxmiber ^ . Equation (69) is used if the final cycle map N was not synthesized. In step 232 
directional information corresponding to 'p is transforaied into three-dimensional coordinates for 
each point of interest Po using, for example, Eqs. (47)-(49) or Eq. (59). 

[0162] It should be appreciated that progressive fiinge division and the combination of 

progressive fiinge division with fiinge synthesis are powerfiil general techniques for resolving 
counting ambiguities and refining counting estimates and that these techniques have applications 
beyond those described above. In particular, these techniques are not limited to the above 
methods and apparatus for acquiring data or even to the analysis of fiinges. 

[0163] With reference to Fig. 16, a variation is now described that utiUzes the 

progressive-fiinge-division and fiinge-synthesis techniques on unresolved fiinges, or image 
speckle. In this figure, laser beam 44 is incident on beam splitter 48 which splits beam 44 into 
beams 46 and 50. Beam 50 reflects off beam splitter 48 and illuminates a region of interest on 
object 10. Beam 46 passes through beam splitter 48 and illuminates reference surface 102. Light 
reflected off of object 10 and reference surface 102 is recombined by beam splitter 48 and passes 
through lens 24 which produces a speckled image on detector array 22. Phase shifting is 
accomplished by inserting a phase-shifting element (not shown) in one of the separated beam 
paths. The frequency of laser beam 44 is variable or can be switched between two or more 
discrete values. Path length differences due to a difference in height z between object 10 and 
virtual image 102' of reference surface 102 produce a modulation of speckle intensity that 
changes with laser firequency. The phase of the modulation at a particular laser frequency can be 
determined accurately through phase-shifting measurements. 

[0164] There are no fiinges produced on the surface of object 10 in Fig. 16 because the 

object 10 is illuminated by only one laser beam 50. Interference occurs at the detector 22 
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because of the coherent superposition of Ught from the two beams 46 and 50. In accordance with 
Eqs. (l)-(5), varying the laser frequency causes the individual speckle lobes in the image of the 
object 10 to oscillate in intensity at a rate that is proportional to the path length difference s. For 
simpHcity, assume that the reference surface 102 is flat and lies in the x-y plane. Let z represent 
the height above the virtual reference surface 102' of the point Po on object 10 to be measured. 
Nominally, s =2z (due to round-trip propagation) and z can be calculated from N by 

if N is obtained through fringe synthesis, and by 



z = ^-^ (72) 
2 V. 



otherwise. 



[0165] The procedure for producing three-dimensional images using Fig. 16 is the same 

as that outlined in Fig. 15, with the exception that in step 232 , Eq. (71), or Eq. (72), is used to 
convert N to z. The approach depicted in Fig. 16 does not rely upon triangulation. 
Consequently, it has the advantage that shadowing effects between source Pi and P2 and receiver 
22 and 24 can be eliminated. This approach is particularly usefiil for shorter-range applications, 
including three-dimensional microscopy, where the separated path lengths remain small. 

[0166] The above embodiments are based on interference using coherent illumination. 

Those embodiments that use a lens 24 to form an image of the object 10 on a detector 22 have 
the potential disadvantage that the resulting image is speckled and that speckle can limit range 
resolution. It should be noted that detectors 23 placed directly on the surface of the object do not 
suffer from this limitation. Speckle occurs for imaging with coherent illumination because the 
finite point-spread ftmction of the lens 24 allows Ught scattered from regions surrounding Po on 
the surface of object 10 to interfere at point Pi at detector 22. The average speckle size at the 
detector plane is determined by the size of the point-spread ftmction of the lens. Speckle can 
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cause a measurement of p to be drawn to one side or the other by scattering regions surrounding 
Po whose contributions can phase up constructively at the image point Pi. 

[0167] A number of techniques are now presented for mitigating the effects of speckle on 

a measurement. In practice, speckle effects tend to average out if many speckles fall within the 
area corresponding to a single pixel of the detector array 22. One way of satisfying this condition 
is to require that the interference fringes are resolved by lens 24 and that the point-spread 
function of lens 24 is small enough to cause many speckles to fall on each element of detector 
array 22. This approach works well in situations that do not require the ultimate lateral 
resolution. Generally, lateral resolution is traded off for speckle averaging. 

[0168] In situations requiring high lateral resolution, it can be necessary to more closely 

match speckle size and pixel size. In general it is desirable to reduce the magnitude of the side 
lobes of the point-spread function of the lens (i.e., to apodize). Doing so reduces the ability of a 
scattering region that would otherwise be contained within a strong side lobe of the point-spread 
function to bias the measurement of p . Apodization can result in a slightly degraded lateral 
resolution. In some instances, it can be advantageous to trade off this lateral resolution for 
improved range resolution. 

[0169] An additional technique for mitigating speckle effects is to minimize the variation 

of ~p over the region of the surface of object 10 corresponding to an individual pixel 22. For 
example, increasing the pixel density while maintaining a constant ratio of speckles per pixel 22 
decreases the change in p from one side of the region corresponding to the pixel 22 to the other 
and improves both range resolution and lateral resolution. 

[0170] For a given optical configuration and pixel density, the variation of p over 

individual pixels 22 can be minimized by making an appropriate choice of the position and 
orientation of the source head. Variation of p is minimized, for example, if the surface element 
corresponding to a pixel 22 lies along an equiphase surface of the three-dimensional fringe 
pattem. In other words, speckle effects are minimized if the surface being measured is parallel 
(or nearly parallel) to a surface of constant phase for the illuminating fringe pattem. This 



-50- 



condition can be achieved by using grazing incidence with the vector representing the orientation 
of the two source points being parallel to the surface normal Except for specialized cases, it is 
impossible to satisfy this condition for every point Po on the surface. The condition can be 
approximated, however, for many objects of practical importance, including, for example, 
aerodynamic surfaces such as auto-body and airframe panels. 

[0171] Many objects have a broad range of surface slopes and orientations, making it 

impossible to satisfy the above condition for minimizing p over the entire surface using a single 
source head. In these situations, this condition can be better satisfied through the use of multiple 
source heads that are positioned such that one of them is responsive to each major region of the 
object 10, or to those regions of greatest interest. It is assumed that some form of multiplexing 
takes place to distinguish which information comes from each source head. If more than one 
source head illuminates the same region, the resulting measurements can be weighted according 
to the level of speckle noise expected in each measurement for that region. Multiple source 
heads can also help to overcome shadowing limitations. Information from multiple source heads 
can be readily combined for a given detector 22 that views a region of overlap. Registration is 
automatic because the angular information corresponding to each pixel is the same for each 
source head. 

[0172] In some instances, it is not possible to illuminate certain regions of an object 10 in 

a way that satisfies the condition for minimizing the variation of p over a pixel. An additional 
technique for mitigating speckle effects that is not affected by this limitation is based on speckle 
averaging. Speckle averaging can be achieved by making additional measurements with slightly 
different parameters (e.g., variations in the optical configurations). Alternatively, the parameters 
can be varied during the exposure. Averaging effectively reduces speckle effects if the set of 
parameters being varied changes enough to produce different statistical reaUzations of the 
speckle pattem. 

[0173] In one embodiment, the mid-position Pm between sources is the adjusted 

parameter. If the direction of displacement of Pm is out of the plane of Figs. 1 and 12, then 
variations of object illumination will be minimized. Variations of Pm are explicitly accounted for 
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in Eqs. (47)-{56) and (59). Reduction of speckle effects can be achieved, for example, by 
averaging the individual measurements of (x, y, z) coordinates for different values of Pm- 

[0174] In another embodiment, the source separation a and the laser frequency v vary 

simultaneously during the exposure time such that the fringe pattem remains stationary. The 
condition for the fringes to remain stationary is seen by inspection of Eq (57) to be that the ratio 
alX (or the product a v) remains constant. To achieve stationary fringes, the motion of the point 
sources Pi and P2 can be slaved to the tuning of the laser 40, or vice versa, by a feedback 
mechanism such as a fringe monitor. The fringe monitor may consist of discrete detectors or a 
linear array of detectors that is incorporated into the source head and measures a fringe pattem 
that is separated from the main beam by a beam splitter. 

[0175] For speckle averaging to be effective, the tuning excursion used is large enough to 

cause a number of decorrelation cycles to occur in the image speckle. Care should be taken to 
minimize the path-length offset So between the separated paths in the source head. Otherwise, 
the phase of the fringe pattem will vary with frequency shift Av by the expression 



A^ = 2;r^ (73) 



Alternatively, the path-length offset can intentionally be set to a nonzero value and the fringe 
pattem "strobed" at 2% increments of the phase Acj) to freeze the phase of the fringe. Different 
phase shifts could then be achieved by varying the strobe timing. 

[0176] It will be appreciated that cycle counting, progressive fringe division, and 

progressive fringe division using fringe synthesis are general techniques for resolving counting 
ambiguities and increasing measurement precision, and that their usefiilness goes far beyond the 
applications described above. Furthermore these techniques can be applied to generalized signals 
(including one-dimensional signals) and are not limited to sinusoidal variations. 

[0177] It should also be appreciated that these techniques apply to general structured- 

light patterns formed by arbitrary means. Examples of means for producing structured light 
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patterns are the interference or diffraction of coherent radiation, the ' Vhite-light" superposition 
of interference or diffraction patterns formed at different wavelengths, and the imaging of 
patterns onto the object with a lens using either coherent or incoherent radiation. As an example 
of the imaging approach, a white-light pattern produced by a digital projector can be made to 
expand and contract to produce accordion motion and to translate sideways to produce phase 
modulation. Although projecting white-light pattems onto the surface in this manner may 
eliminate speckle effects, this approach has certain disadvantages. With an imaged light pattern, 
there is a limit to how small the fringe period may be. It is also necessary to consider the depth 
of field of the projected pattern. It may not be possible to project a fine fringe pattern from an 
oblique angle that is in focus over the entire object. On the other hand, fringes produced by 
interference or diffraction are always "in focus" on the object surface, no matter what the 
distance to the observation point or the size of the fringes. 

[0178] A fiirther embodiment of the invention maintains the advantages of forming 

fringes by interference or diffraction and eliminates speckle effects by forming these fringes in 
broad-band illumination or white light. This embodiment is based on the incoherent 
superposition of interference-fringe intensities formed at different component wavelengths in the 
illumination spectrum. If these fringes have the same size and shape, the interference-fringe 
intensities corresponding to different wavelengths reinforce each other at every point in the 
pattem. 

[0179] By reference to Eq. (57), fringes of the same size and shape are produced by 

fixing the ratio alX. This objective can be achieved simultaneously for each component 
wavelength by using a diffraction grating followed by a lens. For small diffraction angles, the 
angular deviation of the first-order beam diffracted by the grating is proportional to the 
wavelength of that beam. The lens converts this angular deviation into a displacement of the 
source Pi or P2 that is also proportional to thus making it possible to maintain a constant alX 
ratio. In this embodiment, the location of source Pi or P2 is smeared out in wavelength. White- 
light phase shifting can be achieved by translating the grating sideways or through the use of a 
white-light phase shifter, based on, for example, the Pancharatnam phase. 
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[0180] 



There are many possible optical arrangements for implementing a source head 



based on these principles. In one arrangement, Bragg gratings are placed in beams 46 and 50 in 
Fig. 10. The gratings (not shown) and prism 92 are positioned such that the two focal spots 
corresponding to the undiffracted beams overlap and the first-order beams are deviated towards 
the base of prism 92. This arrangement causes the source separation a corresponding to the 
diffracted beams to be proportional to the optical wavelength X. White-light phase shifting, is 
achieved by translating one of the Bragg gratings laterally. Variation of fringe size is achieved 
by changing the period of the grating. In one embodiment, the grating period is changed 
electronically using a Bragg cell, or acousto-optic modulator. In this embodiment, both Bragg 
cells are driven at the same frequency and phase shifting is accompUshed by shifting the phase of 
one drive signal with respect to the other. (The fringes are stationary because each beam suffers 
the same Doppler shift.) In another embodiment, Bragg gratings of different periods are moved 
into place sequentially. In yet another embodiment, a single diffractive element is placed in 
beam 44. A separate white-light phase shifter (not shown) is placed in one of the beams 46 or 
50. In some embodiments of a white-light phase shifter known to those skilled in the art, beam 
splitter 48 is a polarizing beam splitter. 

[0181] Referring now to Fig. 17, therein is shown an embodiment of a broadband or 

white-light interference fringe projector. Source 250 generates a substantially coUimated beam 
of radiation 252 that is directed to diffraction grating 254 at substantially normal incidence. 
Diffraction grating 254 has grating period D. Input beam 252 is represented by three constituent 
wavelength components Xu and X3 for illustration. Beam 252, in actuality, can have arbitrary 
spectral composition. Diffraction grating 254 splits beam 252 into multiple diffracted beams 
whose diffraction orders can be represented by the integer m. For illustration purposes, only rays 
along the perimeter of beam 252 are shown. These diffracted beams propagate at angles 9m with 
respect to the optical axis 258 according to the grating equation for normal incidence, which is 
given by: 




(74) 
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[0182] In one embodiment, diffraction grating 254 is designed to maximize and equalize 

the diffracted efficiency for the diffraction order m = +1 diffracted beam 260 and the diffraction 
order m = -l diffracted beam 262. In other embodiments, diffraction grating 254 is designed to 
maximize and equalize the diffracted efficiency for any set of positive and negative beams of 
equal order and to minimize the energy diffracted into all other orders. Any residual 
undiffracted (m = 0) beam 264, will pass undeviated through diffraction grating 254 and is 
focused by lens 266 onto focal spot 268. 

[0183] The spectral components X,l, A,2, and X3 of focal spot 268 substantially overlap. 

Focal spot 268, in one embodiment, may be substantially blocked by the central obstruction 270 
of optional double slit 272. The different spectral components X\, X2, and of diffracted beams 
260 and 262 are focused by lens 266 onto spectral regions 274 and 276. The distance a{X) 
between the focal spot within spectral region 274 and the focal spot within spectral region 276 
corresponding to a given wavelength X is substantially proportional to the wavelength X. 
Aperture stop 278 of lens 266, in one embodiment, can be used to block undesired higher-order 
diffracted beams. Any undesired residual diffracted orders that pass through lens 266 can be 
blocked, in another embodiment, by the opaque regions of optional double slit 272. Radiation 
from the two spectral regions 274 and 276 expands and overlaps as it propagates and forms a 
broadband interference-fringe pattern 280. Fringe pattern 280 has representative fiinge period d 
at representative distance R from double slit 272. 

[0184] Prior equations such as Eq. (57) which are based on a specific source separation a 

and a specific wavelength X, can be applied to first order to the broadband interference fiinge 
projector illustrated in Fig. 17 by making the substitution 



^ = ^ (75) 

where /is the focal length of lens 266 and D is the period of diffraction grating 254. 

[0185] It should be noted that the fringe-generation scheme depicted in Fig. 17 can also 

produce fringes using narrow-band or laser illumination. One advantage of using diffraction 
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grating 254 followed by lens 266 for narrow-band illumination is that fringe period d is 
insensitive to wavelength so that frequency drifts of the source do not substantially degrade 
measurements. For example, although laser diodes are relatively inexpensive and readily 
available sources, they have a temperature-dependent operating wavelength. However, since this 
technique is insensitive to temperature-dependent wavelength shifts, laser diodes can be used 
without measurement degradation. 

[0186] In one embodiment, diffraction grating 254 is a thin phase grating having a 

square-wave phase profile whose relative phase delay altemates between 0° and 180° for a 
representative wavelength, Xj, with a 50% duty cycle. Grating 254 is relatively efficient, 
diffracting approximately 40.5% of the available energy into each of the m = -1 and m = +1 
diffracted orders, and nominally 0% into the m = 0 and other even diffracted orders. The relative 
phase delay of grating 254 is a fimction of wavelength, causing the energy in the undiffracted 
beam 264 at order m = 0 to increase for wavelengths that differ from the representative 
wavelength X2. In one embodiment, grating 254 is made relatively thin which reduces the 
wavelength dependence of the grating so that it performs well over a wide frequency spectrum. 

[0187] Phase shifting the resulting broadband (or narrow-brand) interference-fiinge 

pattern 280 is achieved by simply translating diffraction grating 254 in the direction 282 shown 
in Fig. 17. White-Hght or broadband phase shifting is realized because a translation of diffraction 
grating 254 by a given fraction of the grating period D shifts each spectral component of fringe 
pattem 280 by the same fraction of the fiinge period d. For example, a translation of grating 252 
by Z)/4, or one-quarter cycle, also shifts the interference-fringe pattem 280 by one-quarter cycle, 
or 90°. 

[0188] Accordion motion (or variation of fringe size) of interference-fiinge pattem 280 

can be achieved in a number of ways. In one embodiment, for small diffracted angles 0m, 
doubling the period D of grating 254 halves the magnitude of Om of beams 260 and 262 (this is a 
small-angle approximation of Eq. (74)), which in turn doubles the period d of fiinge pattem 280. 
In another embodiment, decreasing the focal length / of lens 266 can increase the period d of 



-56- 



fringe pattern 280. The period d of fringe pattern 280 is related to the period D of diffraction 
grating 254, the focal length / of lens 266, and the propagation distance R to first order by 

rf=*5 (76) 

2/ 

[0189] hi order to produce, a continuous variation in d of fringe pattern 280 or to produce 

an arbitrary set oid values lying along a continuum, it is necessary to be able to continuously 
vary one or more of the parameters in Eq. (76). hi one embodiment, lens 266 may be a zoom 
lens with variable focal length / hi another embodiment, the period D of grating 254 may be 
varied by stretching the grating with an applied force, hi yet another embodiment, the period D 
of grating 254 may be varied electronically through, as an example, the variable grating mode 
(VGM) effect in liquid crystal devices, hi yet another embodiment, the period of grating 254 
may vary along the axis orthogonal to the grating lines so that sliding grating 254 sideways in the 
direction perpendicular to 282 can control the period D of grating 254. By using fringe synthesis, 
previously described, it is possible to achieve large variations in effective period d with small 
variations in /, D, or R. For example, fringe-synthesis Eqs. (64) and (67) can be modified as 
follows: 
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(77) 



wrapped 



where Nwrapped is now a ftinction of the ratio fD, 

[0190] It is typically adequate to have a predetermined discrete set of values for / and/or 

D and to resolve fringe-number ambiguities through progressive fringe division as previously 
discussed. In one embodiment, progressive fringe division is accomplished by sequentially 
moving a discrete set of diffraction gratings 254 or lenses 266 into respective positions in Fig. 
17. In alternative embodiments, grating 254 or lens 266 is switchable and assumes discrete 
values of D and / Some examples of switchable gratmgs include liquid-crystal spatial Ught 
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modulators, patterned alignment liquid-crystal spatial light modulators, electrically switchable 
diffractive optical elements, cascaded electrically switchable dif&active optical elements, and 
microelectromechanical systems (MEMS). 



[0191] In one embodiment, illustrated in Fig. 18, diffraction grating 254 is a hquid- 

crystal spatial light modulator 290 consisting of a linear array of long parallel phase retarders 
292, or pixels, whose phase retardance can be controlled to produce a cyclical phase-retardance 
pattem 294. Li this manner, a thin phase grating can be produced and actuated electrically. For 
example, it is possible to produce a square-wave phase profile whose nominal relative phase 
delay altemates between O*" and 180° with a 50% duty cycle. The period D and offset 282 of the 
resulting diffraction grating 254 are variable in discrete steps with high repeatability. Thus, both 
the fringe spacing d and the phase shift of the resulting fringe pattem 280 is controllable 
electrically in discrete repeatable steps at high speeds. 

[0192] In Fig. 18, Wp is the width of pixel 292, and Wbiock is the number of pixels 292 that 

are repeated at the same phase retardance level to produce one half of a 50%-duty-cycle thin 
phase grating. As an example, Wbiock = 2 in Fig. 18. The grating period D is given by 

D = 2«,.^,>v^ (78) 

and can be varied by changing the value of Wbiock- The smallest phase shift Acj) of diffraction 
grating 254, and hence of interference-fringe pattem 280, that can be produced for a given value 
ofwbiock is 



A^=i^ (79) 

'^block 

Thus, for example, phase shifts of ±90**, ±180^, and ±270** can be produced for an «biock value of 2 
and phase shifts of ±60^ ±120^ ±180^ ±240^ and ±300** for an Wbiock value of 3. A generaUzed 
from of Eq. (6) that is valid for three arbitrary phase shifts (|)2, and (t)3 can be expressed as: 
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^ ^(cos^2 " 003^3)+ /2(cos<;»3 -cos(»^)+/3(cos<^^ -cos^^) ^^^^ 
/i(sin^3 -sin^2)+A(siii^^i "Sin(^3)+/3(sin^2 -sin^,) 

where lu h, and h are signal intensities corresponding to phase shifts (1)2, and <|)3. Thus, the 
phase shift values can be selected differently for each value of Wbiock according to the availability 
of phase-shift values. 

[0193] Fringe synthesis is an effective means of obtaining additional values of D beyond 

those achievable by varying the integer Wbiock in Eq. (78). The available values of the effective 
grating period Aff that can be synthesized from two measurements with different values of Wbiock 
are given by 



^e.= 1 1 (81) 



'^blockl '^block2 



where Wbiocki < «biock2. Large values of Deff can be obtained using Eq. (81) by making Wbiocki and 
Wbiock2 successive integers. As an illustration of a possible measurement strategy, assume that the 
pixel width W;, is 10 |im and that a sequence of three measurements is taking with Wbiocko = 2, 
Wbiocki = 19, and «biock2 = 20. This sequence can be used to produce a first grating period D of 40 
\\m using «biock = 2 in Eq. (78) and a second grating period of 400 \xm using Wbiock = 20 in Eq. 
(78). A third, and much larger grating period, of 4.2 mm can be synthesized by using Wbiocki = 19, 
and «biock2 = 20 in Eq. (81). The ratio between successive grating periods in this sequence is 
approximately ten, which is a reasonable sequence for progressive fringe division. 

[0194] Referring back to Fig. 1 7, wide variations of nbiock produce corresponding wide 

variations in the spread and location of spectral regions 274 and 276. In one embodiment, a more 
sophisticated version of double slit 272, such as beam block 296 shown in Fig. 19, may then be 
used to eliminate unwanted residual diffracted orders. Beam block 296 comprises two nested 
double slits and includes central obstruction 268. The outer double slit 275 passes spectral 
regions 274 and 276 for the smallest value «biocko used in the measurement sequence in the above 
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example, while the inner double slit 277 passes the desired spectral regions 274' and 276' for the 
larger values Wbiocki and Wbiock2. 

[0195] In yet another embodiment illustrated in Fig. 20, accordion motion is achieved by 

placing lens 300 on optical axis 258 such that images 302 and 304 of spectral regions 274 and 
276 are formed at image plane 306. Note that in another embodiment, images 302 and 304 may 
be virtual images (not shown), if the focal length of lens 300 is negative. Accordion motion 
results in one embodiment, by varying the magnification of images 302 and 304 by moving lens 
300 along optical axis 258. Or, in another embodiment, accordion motion results by changing 
the focal length of lens 300. Li yet another embodiment, lens 300 is a zoom lens, which varies 
the focal length without changing the location of image plane 306. In still an altemative 
embodiment, a discrete set of lenses of different focal lengths can be placed in suitable positions 
sequentially. The use of lens 300 produces accordion motion without changing any of the 
components within Fig. 17, including double slit 272. 

[0196] Yet another embodiment of a broadband interference-fringe projector is shown in 

Fig. 21. This embodiment includes broadband soxirce 250, an acousto-optic modulator (AOM) 
306, and a lens 266. AOM 306 is driven with a compound signal 308 generated by multiplier 
310 which multiplies a sinusoid 312 having a variable frequency /„ by a sinusoid 314 having a 
fixed frequency fc. Compound signal 308 produces two traveling acoustic waves of different 
periods in AOM 306 that in tum generate two angularly separated beams 260 and 262. The 
different spectral components of diffracted beams 260 and 262 are focused by lens 266 into 
spectral regions 274 and 276, where the separation aQC) between the focal spot within spectral 
region 274 and the focal spot within spectral region 276 corresponding to a given wavelength X 
is substantially proportional to the wavelength X. In addition, the midpoint between like spectral 
components of spectral regions 274 and 276 is substantially localized at position 268 for each 
spectral component. By adjusting^m the separation fl(A.) for each wavelength component X is 
adjusted about point of symmetry 268 from zero separation to the maximum separation permitted 
by the bandwidth of AOM 306. Radiation from the two spectral regions 274 and 276 expands 
and overlaps as it propagates to form a broadband interference-fringe pattem 280 of 
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representative period d at representative distance R. Varying/,, causes fringe pattern 280 to 
expand or to contract. 

[0197] Since the pressure variations across AOM 306 are traveling, beams 274 and 276 

are Doppler shifted with respect to one another. This causes fringe pattern 280 to travel across 
the region of overlap at a speed proportional to the difference frequency between the two 
frequency components of compound signal 308. hiterference pattern 280 can be effectively 
frozen at a particular spatial phase by amplitude modulating source 250 at the difference 
frequency 2fm of the two beams 260 and 262. The spatial phase of the fringes 280 can be 
precisely controlled by changing the phase (e.g. changing the delay) of the laser drive signal 316 
with respect to the AOM control signal 308. Alternatively, drive signal 316 can be used to 
control an external amplitude modulator placed in the optical path (not shown) or to modulate the 
sensitivity of detector 22 used to observe the fringe pattern. 

[0198] In the above embodiments of a broadband interference-fringe projector, beam 252 

is substantially coUimated so that diffractive element 254 (Fig. 17) or 306 (Fig. 21) diffracts each 
spectral component of a given diffraction order m into the same diffracted angle 9^. Other 
variations of this configuration are also possible, and remain within the spirit and scope of this 
invention. For example, in alternate embodiments, beam 252 may be either diverging or 
converging. In yet another embodiment, beam 252 may enter diffractive element 254 at oblique 
incidence. In still another embodiment, the order of the refractive and diffractive elements may 
be reversed. 

[0199] Many alternative Ught sources 250 may be used, depending on the degree of 

temporal and spatial coherence desired. Low temporal coherence (broad spectral content) 
reduces or eliminates speckle effects. Whereas, high spatial coherence allows for tighter focal 
spots to be formed for each spectral component in spectral regions 274 and 276. Generally, tight 
focal spots improve the visibility of fringes and increase the depth of field of fiinge pattem 280. 
Note that the depth of field is infinite for spatially coherent sources. Additionally, tight focal 
spots also make it easier to separate spectral regions corresponding to different diffracted orders 
m so that undesired spectral regions can be blocked with representative masks 272 or 296. 
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[0200] One method of achieving high spatial coherence is to make light source 250 laser 

based. In one embodiment, source 250 is a laser diode with collimating optics. Laser diodes 
with spectral widths of approximately 3nm are commonly available. In another embodiment, the 
spectral width can be broadened by coupling together an array of laser diodes operating at 
different wavelengths into a single beam. In still another embodiment, a rapidly scanning tunable 
laser can be used. Another approach to broadening the spectrum of a laser is to pump a standard 
single-mode silica fiber with a passively Q-switched microchip laser. A variety of nonlinear 
effects in the fiber produce an extremely broadband spectrum at the output end of the fiber. In 
yet another approach, broadband illumination is achieved with a modelocked titanium:sapphire 
laser. 

[0201] There are many other approaches for generating beam 252 that rely on more 

conventional sources. Lower (but adequate) spatial coherence can be achieved, for example, by 
coupling an arc lamp into an optical fiber. The output end of the fiber acts like an extended 
source. If the output of the fiber is coUimated using a lens, the longer the focal length of the lens, 
the smaller the variation in the angle of propagation in beam 252 and the tighter the focal spots 
that can be formed within spectral regions 274 and 276. When using an extended source, spatial 
incoherence effects are minimized by placing lens 266 at a distance along optical axis 258 such 
that the surface of diffraction grating 254 or the active area of AOM 306 would be imaged onto a 
plane on or near the surface of object 10 being illuminated. This configuration minimizes spatial 
incoherence effects because the rays that reach a given area of fiinge pattem 280 to interfere all 
originate from the same region of the dififractive element where the input beam is spht into 
different diffracted orders. 

[0202] The above embodiments for projecting fiinge patterns (whether based on narrow- 

band or broad-band illumination) vary in the accuracy and repeatability to which the phase shifts 
(f>\, ^, and (fh, can be set. For instance, phase shifting by moving an optical component such as 
mirror 94 or 96 in Fig. 10 or grating 254 in Fig. 17 with a piezoelectric transducer suffers from 
inconsistencies due to hysteresis. Variations of these phase shifts from their assumed values 
cause errors in the calculation of (j) in Eq. (6) or Eq. (80). These errors in ^ are periodic in nature 
and produce a rippled appearance in the reconstructed surface profile. Likewise, a rippled 
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appearance is caused by imperfections in the illumination sequence. These imperfections may be 
caused, for example, by changes in the exposure level or bias level during acquisition of signal 
intensities /i, and /a. Exposure-level variations may occur, for instance, from fluctuations in 
source intensity or from differences in exposure time due to shutter imperfections. 

[0203] The ripple pattem occurs on the reconstructed surface with the same orientation as 

the projected fringe pattem on the physical surface 10. The dominant period of the ripple pattem 
depends on the error source. Phase-shift errors, for example, produce two ripple oscillations per 
fringe period. Ripple can be the dominant error source in a surface profile measurement, 
particularly if speckle has been eliminated through broad-band illxmiination, or if speckle effects 
from narrow-band illumination have been suppressed through one of the mitigation techniques 
described above. 

[0204] A versatile and robust technique, referred to as histogram-based derippling, is now 

introduced that corrects for ripple error. In this technique, the data is analyzed to determine actual 
phase offsets and exposure-level variations so that ^can be calculated using corrected values. In 
order to elucidate this technique, Eq. (80) is first rewritten in a form where assxmied values, as 
well as quantities calculated from assumed values, are differentiated from actual values: 

^ /^|cos^2 -cos^3y /;^os4 -cos^^^/g^os^^- cos^3^ 

(jf = tan z'* z*^ '*\ / ^ *\* (^^) 

/|fsin^3 -sin^2 Vh/jfein^^i -sin^j W/jfon^z -sin^ j 

In Eq. (82), ^ , ^2 > and ^3 are the assumed phase-shift values and ^ is the calculated phase using 
these assumed values. 

[0205] Variations due to imperfections of the illumination sequence may be incorporated 

into Eq. (82) by writing the intensity in the form 

. (83) 

where the subscript / takes one of the values 1, 2, or 3, and represents a particular member of the 
illumination sequence. The parameters Ai and Bi allow for variations in the exposure level and 
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bias level of the illumination sequence, respectively. Without loss of generality, the nominal 
peak-to-valley amplitude swing in Eq. (83) can be assumed to be unit valued. 

[0206] If the values of Ai and 5, are known, or are estimated as ^, and 5, , imperfections 

in the illumination sequence can be corrected by using the calculated intensity 

/,=^ (84) 
A 

to replace the measured intensity/,- in Eq. (82). In addition, phase-step errors can be eliminated if 
the actual values ^, and ^ can be determined and used in Eq. (82) rather than the assumed 
values ^ , ^2 > ^d ^3 . One approach to minimizing ripple is to vary the parameters Ai, Bi and (jn 
such that a metric indicating the degree of ripple in ^ is minimized. Besides being 
computationally intensive to implement and requiring a multi-dimensional search algorithm to 
minimize the metric, it is difficult to find a suitable metric for minimization without first 
unwrapping the phase function ^ over a representative subset of the surface. One example of a 
metric that could be used, once unwrapping has taken place, is the standard deviation of the 

A A 

difference between^ and a smoothed representation of ((>. This metric tends to decrease as the 
contributions due to ripple decrease. 

[0207] Histogram-based derippling is superior to the above derippling approach because 

it does not require phase unwrapping and produces corrected parameters directly. Reference is 
now made to Fig. 22 for a continued introduction to this technique. Figure 22 contains plots of 
the error (^>-^)/(27r) as a function of the calculated quantity ^/(27i) for different values of the 
parameters Au Bi, and ^. This error is equivalent to the ripple-induced component of the fiinge- 
number error AA^ and limits range resolution as described in Eq. (32). In these curves, the 
assumed phase-step values are ^ =-120°, ^2 = 0°, and ^3 = 120°. Both ^and ^ are considered 
to be wrapped quantities that range between -71 and + tx so that ^{2n) and ^/(2n) range between 
-1/2 and +1/2. Curve (a) shows the effects of amplitude variations A\ = 0.9, = 1, ^3 = 1.2, 
Curve (b) corresponds to bias variations B\ = 0.2, ^2 = 0, ^3 = 0.1, and Curve (c) displays the 
effects of having actual phase shifts ^ = -105°, <fc = 0°, and ^ = 100° that differ fi-om the 
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assumed values. Finally, Fig. 22a is a plot of the slope of the curves in Fig. 22, given by the 
derivative of ^ with respect to ^ and denoted as ^(^). This slope oscillates about unity and 
converges to unity at all points as the assumed parameters approach their true values. Note that 
relatively large slope variations are produced for relatively small errors in AN, making slope a 
sensitive indicator of error in the assumed values of Ai, Bi, and ^. 

[0208] Clearly, Fig. 22 and Fig. 22a contain a great deal of information about Ai, Bi, and 

<f>i. For example, these plots illustrate that bias variations introduce one ripple per fringe and that 
phase-step errors introduce two ripples per fringe. Amplitude variations produce a mixtwe of one 
and two ripples per fringe. If it were possible to produce the curves graphed in Fig. 22 or Fig. 22a 
from analysis of the data, the dependence of the governing equations for these curves on Ai, Bi 
and ^ could be used to recover information about these parameters. Then the data could be 
derippled by calculating fusing these corrected parameters. The difficulty with this approach is 
that these graphs presuppose knowledge of the true values of <j>, 

[0209] The histogram approach overcomes this difficulty. In this approach, a histogram is 

first made representing the relative occurrence of each value of ^ in the measurement (or some 
subset of the measiwement). The histogram is created by dividing the range of possible values of 
^ into a finite number of bins and counting the number of occurrences of ^ that fall within each 
bin. The histogram of ^ provides a numerical approximation to the probabiUty density fiinction 
of (j), written as The two probability density functions P(^) and P{(^ can be related 
through the well-known relation 

P{^)= P{^)\<l>'{^)\ , (85) 

which assumes that ^ is a single-valued function of ^ . The derivative in Eq. (85) is what is 
graphed in Fig. 22a. The probability distribution function P(^) is available from the data through 
calculation of the histogram of ^. Although P(^ in Eq. (85) is unknown, it is reasonable to 
assume for most cases of interest that is xmiformly distributed over its range of potential values. 
This assumption becomes more precise as the number of fiinges covering the sampled surface 
area increases. In practice, the assumption that P{(^ is uniformly distributed has proven to be 
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adequate for even a relatively small number of fringes. Other forms of however, can be 
used for histogram-based derippling where the functional form of is refined based on 
additional information about the surface. 



[0210] 



A 

Given P(<zJ), Eq. (85) can be used to relate the quantity P(^) (obtained through 



calculating a histogram of the measured values ^) to the absolute value of the derivative <*'(^) • 



Actual values of Ai, Bi and (fn can then be determined based on the functional form of this 
derivative. To illustrate how this may be done, the functional form of <*'(^)is now analyzed in 
greater detail 



Bi = 0), but only phase step errors. For this case, it can be shown that <l>'{^) takes the functional 
form 



where K, Kd, and Ksz are constants whose values depend on the assumed phase-step values (/>i 
and the actual phase-step values ^ . The value of AT (which is close to unity) can be determined 
from the other two constants Kc2 and Ksi, leaving only two degrees of freedom in Eq. (86). The 
constants Kc2 and Ks2 may be obtained from the shape of the histogram of ^ by, for example, 
taking the reciprocal of this histogram and calculating the Fourier coefficients that correspond to 
the cos(2^) and sin(2^) terms of a Fourier series expansion. The Fourier coefficients for the 

A A 

cos(2^) and sin(2^) terms are then normalized by the DC term of the Fourier series to yield Kci 
and Ks2, respectively. 

[0212] Because Eq. (86) only has two degrees of freedom, it is not possible to fully 

recover the values of the three phase shifts ^, ^, and ^ from this equation. Instead, what can be 
recovered are the relative phase shifts and two of which determine the third. 

Knowledge of these relative phase shifts is sufficient to completely eliminate ripple. There is. 



[0211] 



First, consider the case where there are no amplitude or bias variations (Ai = 1 and 




(86) 
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A 

however, a small offset error in the improved calculation of ^ that remains. This offset error is 
due to lack of knowledge of the true value of ^. 



[0213] For simplicity, let the assumed phase shifts be equal, of magnitude (j^tep, and 

centered about zero. These phase shifts can then be expressed as^i = -^^/e^, (j>2 = 0, and ^3= ^^^ep. 
Also assume that ^ = 0 by default. The values of ^ and ^ can then be determined from Ka and 
Ks2 by the relations 



^ =-cos 



(1 + 2^2 +^-=2 cos^„<5p + A;^ sin^„^) 



(87) 



and 



= COS 



^ , -hi^^^, + (u /r,, >os - (1 -h ir,,X2 sin ^ 
(1 + i:,2 Xl + cos^,,^ - sin ^,,^) 



. (88) 



For the case previously assumed in Eq. (6) where (l>step = 90°, Eqs. (87) and (88) reduce to 



^ =-cos 



( 2 2^ 

K^cl K2 + + K2 + f^clf^sl 



(89) 



and 



^3 = cos ' 



/ 2 2 ^ 

f^c2 + — K^jK. 



0+^.2X1-^.0 



(90) 



[0214] As a further example of the application of histogram-based derippling, the 

technique is applied to a combination of phase-step errors and amplitude errors. Although the 

A 

expression for ^(^) for this case is much more complicated than Eq. (86), it can be 
approximated for small values of land ^— ^, as 



1 



(91) 
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[0215] Once again, the coefficients Kcu Ksh Ka, and Ksi are obtainable through 

Fourier analysis, being the corresponding Fourier coefficients of a Fourier-series expansion of the 
reciprocal of the histogram used to approximate ^'(^) . 

[0216] Setting ^2 = <fe = 0 and = 1, and allowing unequal step sizes for the assumed 

values of ^ and ^3, the newly approximated values of ^, ^, i4i, and can be determined from 
Kcu Kc2y and Ks2 through the relations 

^ = ^ + a;, sin^ +i:„(l -cos^,y A;2^in^3 -sin(^, +^3)^ A:,2(cos^3 -cos(^, + ^3)^ , (92) 

A A ^ ^ \ / AAX ^ ^ AAV 

^3 = ^^3 + a;, sin ^3 + ^1 - cos^3 j+ K,^ ^sin ^ - sin(^, + ^3 ) K,^ ^os ^ - cos(^i + ^3 ) J . (93) 

A, =l + iCi(l- cos^)-^:,, sin^. , (94) 
and 

^3 = l+^ci(l-cos^3)-/:,,sin<^3 . (95) 

[0217] Because Eqs. (92)-(95) are based on the small-error approximation given by Eq. 

(91), the calculated values of ^, ^, Au and A^ are improved, but not exact. If desired, higher 
accuracy can be obtained by iterating the procedure. To do so, the values of ^ and ^ calculated 
from Eqs. (92) and (93) become the new assumed values ^ and ^3 in Eq. (82), and the values of 
A\ and A^ obtained from Eqs. (94) and (95) become the new assumed values and in Eq. 
(84). This procedure converges rapidly, providing many significant decimal places of precision in 
only a few iterations. 

[0218] As an additional benefit, bias variations Bi are handled automatically within the 

context of this iterative approach. As the number of iterations increases, the effects of bias 
variations Bi are eliminated by being folded into the other parameters Ai and ^. This approach 
eliminates ripple, leaving a small offset error similar to the error introduced by the lack of precise 
knowledge of ^. 
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[0219] Having described and shown the preferred embodiments of the invention, it will 

now become apparent to one of skill in the art that other embodiments incorporating the concepts 
may be used and that many variations are possible which will still be within the scope and spirit 
of the claimed invention. It is felt, therefore, that these embodiments should not be limited to 
disclosed embodiments but rather should be limited only by the spirit and scope of the following 
claims. 
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